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Abstract 


A  method  for  the  simultaneous  transmission  of  multiple  sonar  signals  is  described 
and  investigated.  The  sonar  signals  are  linearly  combined  as  the  modulating  phcise  com¬ 
ponent  of  a  carrier  signal,  which  avoids  the  transmit  energy  losses  in  peak  power  limited 
systems  caused  by  a  modulating  amplitude.  These  modulating  signals  are  designed  to 
be  spectrally  separable  upon  reception,  so  each  sonar  signal  can  be  simultaneously  pro¬ 
cessed.  Demodulation  and  echo  processing  techniques  are  then  presented,  and  several 
experiments  (both  simulated  and  in-water)  were  conducted  to  investigate  the  detection 
performance  of  the  modulated  signals,  as  compared  to  the  equivalent  unmodulated  sig¬ 
nals.  Results  from  these  experiments  indicate  that  the  modulation  technique  is  a  viable 
method  of  transmitting  multiple  sonar  signals  simultaneously,  while  avoiding  amplitude 
modulation. 
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Chapter  1 


Introduction 


1.1  Motivation 

The  detection  of  a  scatterer  and  the  estimation  of  it’s  parameters  are  some  of 
the  main  goals  of  active  sonar.  Sonar  waveforms  are  transmit  into  the  ocean,  and  their 
echo  signals  are  processed  for  information  pertaining  to  range  and  velocity  characteristics 
of  any  scatterers,  that  the  signals  encounter  [36].  The  ability  to  successfully  estimate 
these  scatterer  parameters  in  different  situations  depends  on  the  sonar  waveforms  used  to 
gather  this  information.  Some  waveforms  have  good  delay  (or  range)  resolution  but  poor 
Doppler  (or  velocity)  resolution,  while  others  exhibit  excellent  Doppler  resolution  but 
inferior  delay  resolution.  Some  waveforms  perform  well  in  the  presence  of  reverberation, 
while  others  are  good  at  distinguishing  between  multiple  scatterers  closely  spaced  in  delay 
and  Doppler.  Since  different  sonar  waveforms  are  better  suited  for  different  situations  and 
estimating  certain  scatterer  parameters,  it  is  beneficial  to  use  multiple  sonar  waveforms 
to  interrogate  a  given  scatterer. 

The  parameters  of  a  scatterer  can  change  over  time.  Although,  the  interrogation 
of  a  scatterer  with  multiple  sonar  signals  over  different  time  intervals  may  have  some 
benefit,  it  is  more  useful  to  simultaneously  transmit  several  waveforms,  because  they 
estimate  the  scatterer  parameters  for  a  specific  event.  Since  transmission  time  is  limited 
in  monostatic  sonar,  the  simultaneous  transmission  of  multiple  sonar  waveforms  also 
maximizes  transmission  time,  which  in  turn  maximizes  the  Doppler  resolution  of  the 
sonar  waveforms  [15]. 

Although,  several  methods  do  exist  to  simultaneously  transmit  multiple  wave¬ 
forms,  some  of  the  most  widely  known  ones  come  from  the  field  of  radio.  They  are 
amplitude  modulation,  which  is  the  modulation  of  a  carrier  waveform’s  amplitude  with 
multiple  waveforms,  and  angle  modulation,  which  is  the  modulation  of  a  carrier  wave¬ 
form’s  phase  with  multiple  waveforms.  While  both  modulation  techniques  provide  viable 
methods  of  simultaneously  transmitting  multiple  sonar  waveforms,  amplitude  modula¬ 
tion  experiences  total  transmit  energy  losses  in  peak  power  limited  systems  [30].  These 
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energy  losses  present  problems  for  amplitude  modulation,  since  scatterer  detection  de¬ 
pends  on  the  transmit  level  [36]. 

Phase  modulation  (a  subset  of  angle  modulation)  provides  a  method  for  trans¬ 
mitting  multiple  sonar  waveforms  simultaneously,  while  maximizing  transmit  time  and 
avoiding  amplitude  modulation.  Phase  modulation  is  the  linear  combination  of  multiple 
waveforms  as  the  modulating  phase  component  of  a  carrier  waveform  [41].  When  the 
modulating  waveforms  are  designed  to  be  spectrally  separable  by  filtering  techniques, 
the  modulating  signals  can  be  individually  processed  upon  reception  and  demodulation. 
The  application  of  this  technique  in  active  sonar  is  referred  to  as  sub-band  modulation, 
which  was  first  proposed  by  L.  H.  Sibul  [30]. 

This  thesis  presents  a  discussion  of  modulation  theory,  as  it  applies  to  sub-band 
modulation.  Difierent  modulation  techniques  are  compared,  and  a  detailed  analysis  of  the 
sub-band  modulation  and  demodulation  processes  is  made.  Furthermore,  the  waveform 
design  considerations  for  both  the  sub-band  waveforms  and  the  overall  modulation  are 
presented,  and  sub-band  modulation  theory  is  implemented  as  a  process  in  Matlab.  Both 
simulated  and  in-water  experiments  are  conducted  to  investigate  the  effects  of  scatterer 
velocity,  sub-band  modulating  signal  type,  signal-to-interference  ratio  and  demodulation 
conditions  on  the  performance  of  sub-band  modulating  signals,  and  to  demonstrate  the 
application  of  sub-band  modulation  in  cictive  sonar. 

1.2  Literature  Review 

The  coupling  of  sonar,  signal  processing  and  modulation  theory  provided  the 
basis  of  sub-band  modulation  theory.  This  is  because  the  concepts  behind  sub-band 
modulation  are  not  new,  only  the  application  is.  Therefore,  a  wide  variety  of  sources 
were  utilized  to  construct  the  theoretical  background  for  this  thesis. 

Numerous  sonar  and  sonar  signal  processing  sources  were  employed  in  this  thesis 
to  provide  a  solid  background  in  sonar  theory  and  the  motivation  behind  sub-band 
modulation.  Some  sources  provided  a  detailed  background  in  underwater  propagation 
and  sonar  theory  [3,  23,  36],  while  others  provided  insight  into  the  benefits  of  different 
sonar  waveforms  [4,  15,  32].  Still  other  sources  provided  am  excellent  insight  into  sonar 
signal  processing  [1,  5,  11,  35,  37,  38]. 

Extensive  resources  exist  on  the  topic  of  modulation  theory.  Since  modulation 
theory  was  exhaustedly  researched  through  out  the  twentieth  century  because  of  it’s 
applications  in  radio,  a  vast  array  of  literature  can  be  found,  that  present  comprehensive 
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overviews  on  both  amplitude  and  angle  modulation  [2,  13,  27,  29,  41].  Other  more 
specific  sources  focus  on  the  effects  of  noise  on  modulation  [9,  19,  28,  40]. 

1.3  Thesis  Outline 

Chapter  2  presents  an  in-depth  background  into  active  sonar  and  the  motivation 
behind  sub-band  modulation.  The  active  sonar  background  provides  insight  into  the 
effects  of  scatterer  parameters  and  underwater  sound  propagation  on  the  echo  signal. 
A  discussion  of  echo  processing  is  then  presented  to  illustrate  the  process  of  estimating 
these  parameters  from  received  echo  signals.  This  is  followed  by  a  comparison  between 
several  sonar  waveforms  and  a  discussion  of  their  advantages  and  disadvantages  in  var¬ 
ious  situations.  Finally,  the  motivation  behind  sub-band  modulation  is  presented,  as  a 
culmination  of  the  ideas  presented  in  the  previous  sections. 

A  comprehensive  discussion  of  modulation  theory  and  it’s  application  to  sub-band 
modulation  is  given  in  chapter  3.  First,  an  overview  of  amplitude  modulation  is  presented 
to  illustrate  it’s  limitations.  An  overview  of  angle  modulation  was  then  discussed,  and 
it’s  subsets  (frequency  and  phase  modulation)  were  compared  to  amplitude  modulation. 
This  provides  a  clear  definition  of  sub-band  modulation,  along  with  an  explanation  of  it’s 
advantages  in  active  sonar  over  other  types  of  modulation.  Next,  specific  characteristics 
of  sub-band  modulation  were  examined  and  discussed.  Chapter  3  ends  with  an  analysis 
of  the  demodulation  of  sub-band  modulating  signals  for  processing.  The  overall  demod¬ 
ulation  process  is  first  presented,  then  the  effects  of  delay  and  Doppler  on  demodulation 
are  analyzed,  and  finally,  demodulation  in  the  presence  of  interference  is  discussed. 

In  chapter  4,  the  experimental  portion  of  this  thesis  is  presented.  First,  a  detailed 
description  of  the  factors  infiuencing  sub-band  waveform  design  are  discussed,  along 
with  the  waveform  design  process  that  takes  them  into  account.  This  is  followed  by  the 
presentation  and  analysis  of  the  sub-band  waveforms  used  in  this  experiment,  and  the 
motivation  behind  each  waveform  combination  is  examined.  Next,  the  implementation  of 
sub-band  modulation  theory  as  a  series  of  Matlab  processes  is  explained.  A  step  by  step 
description  of  the  simulations  in  this  thesis  is  presented  in  the  following  section,  which 
focuses  on  the  setup  of  the  different  experiments  and  their  goals.  A  similar  description 
of  the  in-water  experiment  portion  of  this  thesis  concludes  this  chapter. 

Chapter  5  consists  of  a  discussion  of  the  results  of  the  experiments  described  in 
chapter  4.  First,  the  results  from  each  of  the  simulations  were  incrementally  presented 
and  analyzed.  This  analysis  is  then  used  to  characterize  the  detection  performance  of  the 
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sub-baad  modulated  signals,  as  compared  to  the  performance  of  equivalent  unmodulated 
signals.  The  results  from  the  in-water  experiment  are  then  discussed  and  analyzed. 
These  results  were  then  compared  with  the  simulation  results  to  confirm  the  detection 
performance  of  the  sub-band  modulated  signals  observed  in  the  simulations. 
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Chapter  2 


Sonar  Theory  and  Motivation 


A  thorough  background  in  active  sonar,  sonar  signal  processing  and  the  motivation 
behind  sub-band  modulation  is  warranted  to  better  understand  sub-band  modulation 
and  its  role  in  active  soneir.  First,  active  sonar  theory  is  presented  in  order  to  illustrate 
the  resulting  effects  of  scatterer  parameters  on  the  received  echo  signal.  Second,  the 
echo  processing  techniques  used  to  estimate  these  scatterer  parameters  from  the  echo 
signal  are  discussed.  Third,  several  sonar  waveforms  are  compared  in  order  to  bring 
to  light  the  advantages  and  disadvantages  of  different  sonar  waveforms.  Finally,  eictive 
sonar  theory,  echo  processing  and  the  characteristics  of  different  sonar  waveforms  are 
combined  together  to  motivate  and  begin  to  define  sub-band  modulation. 

2.1  Active  Sonar  Background 

Active  sonar  is  the  interrogation  of  the  ocean  environment  with  acoustic  signals 
to  detect  the  presence  of  possible  scatterers  from  backscatterered  echos  and  to  estimate 
scatterer  parameters,  such  as  velocity  and  range  [3,  15,  36].  A  sonar  waveform  has  the 
form 


fit)  =  [«(*)  exp(jV’(i))]  expOwgt) 

=  /(t)  expOwpt)  t€(0,T),  (2.1) 

where  / (t)  denotes  the  complex  envelope  of  the  waveform  and  contains  all  of  the  phase 
and  amplitude  information.  The  analytic  representation  of  the  real  waveform  is  given  by 

^{fit)}  =  a(i)  cos  )  (2-2) 

where  is  the  frequency  of  the  carrier  term,  a{t)  represents  the  amplitude  of  the 
signal  and  'tpit)  represents  the  phase  [5,  11,  35,  39].  Sonar  waveforms  by  convention  are 
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normalized  to  unit  energy, 


-oc  0 


[30,  35,  38],  so  that  the  transmitted  signal  with  energy  is 


»,(<) 


*!{/(<)) 


^ye^a{t)  cos  +  ipit)  . 


(2.4) 


The  received  signal  is  the  combination  of  the  echo  signal  and  additive  interference, 
consisting  of  white  Gaussian  noise  and  reverberation  measured  at  the  receiver  [1,  11,  35]. 
This  signal  is  modeled  as 


r{t)  =  s{t)  +  n{t), 


(2.5) 


where  s{t)  is  the  echo  signal  and  n{t)  is  the  interference.  The  simplest  backscatter  model 
for  the  returning  echo  signal  is  the  point  reflector  in  a  homogeneous  environment,  where 
the  echo  signaJ  is  an  attenuated,  time  delayed  and  possibly  Doppler  shifted  replica  of 
the  transmitted  waveform  [15,  23,  35].  The  echo  signal  is 

s{t)  =  I  ^')  ^  [^'  ~  ^0] 

0  ^ 


where  [p' is  the  echo  spreading  function  [15,  22).  The  echo  spreading  function 
is  the  random  weighting  of  the  echo  signal,  that  accounts  for  the  scatterer-influenced 
attenuation,  delay  and  Doppler  changes  in  the  echo.  for  a  point  scatterer 

with  dilation  and  delay  r  (commonly  known  as  a  slowly  fluctuating  point  scatterer)  is 

=  (2.7) 


[23,  38],  where  5^  has  a  random  Rayleigh  magnitude  with  uniform  phase  that  remains 
constant  over  each  interrogation  interval  ^5^  =  exp  (jO)J  [15,  35].  (2.5)  can  be 
rewritten  to  incorporate  the  spreading  function  of  the  slowly  fluctuating  point  scatterer 
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using  (2.7)  as 

r{t)  =  j  j  yfj' SJ  (0'  -  <5  (r'  -  /  ^0'  (t  -  /)]  dr  +  n{t) 

=  ^SJi0{t-r))  +  n{t) 

=  ^SJ{p,t)  +  n{t),  (2.8) 

where  /  (p,  f)  is  the  wideband  Doppler  model  of  the  echo  waveform,  and  p  represents  the 
delay  and  Doppler  parameters, n{t)  is  white  Gaussian  noise  with  a  power  spectral 
density  of  N^,  when  noise  is  the  dominant  form  of  interference  [23].  Reverberation 
can  dominate  the  interference,  and  this  possibility  will  be  further  discussed  later  in  this 
section.  The  magnitude  of  the  echo  signal  is  also  affected  by  equipment  and  transmission 
medium  influenced  parameters.  These  attenuation  parameters  can  be  incorporated  into 
the  variable  S  ,  along  with  target  strength. 

3 

The  wideband  Doppler  model  of  the  received  echo  signal  can  be  rewritten  in  the 
form  of  (2.1),  where 

=  yle^Sj^{0{t-T))exp(juj^lpt-Tfj,  (2.9) 

r  is  the  time  delay,  and  0  is  the  dilation  factor,  represented  by 


«  1  +  -^.  (2.10) 

c 

[1,  22,  38].  V  is  the  platform  velocity,  v  is  the  scatterer  velocity,  and  is  the  relative 
line-of-sight  velocity  between  the  platform  and  scatterer,  which  is  positive  when  the 
range  between  the  platform  and  the  scatterer  is  closing  and  negative  when  it  is  opening. 
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The  receiving  and  transmitting  transducers  are  colocated  for  a  monostatic  sonar  [36],  so 
that  the  round  trip  delay  is 


2R 


(2.11) 


where  R  is  the  distance  to  the  scatterer  and  c  is  the  speed  of  sound  [15,  35,  39].  The 
range  and  scatterer  velocity  are  estimated  by  measuring  the  dilation  factor  and  time 
delay  of  the  echo  signal. 

Equation  (2.9)  illustrates  the  compression  or  dilation  of  the  time  scale  of  the 
complex  envelope  and  the  frequency  shift  of  the  carrier  frequency.  It  can  be  rewritten  as 


s{t)  =  yJe^Sj{0{t  -  t))  expiju^it  -  t))  exp(j27r^t)  (2.12) 

by  separating  the  delay  and  dilation  components  of  the  carrier.  The  shift  in  the  carrier 
frequency  caused  by  the  motion  of  the  scatterer  relative  to  the  sonar  platform  is 


and  is  known  as  the  Doppler  frequency  [35,  38]. 

The  signal  (2.12)  may  be  simplified,  provided  the  scatterer  does  not  travel  a 
greater  distance  than  the  range  resolution  during  the  time  of  its  illumination  by  the 
sonar  signal.  This  requires  that 


c 


(2.14) 


or  •C  c  to  unambiguously  resolve  a  scatterer  [25,  32,  35,  38].  Here,  T  is  the  signal 
duration,  and  W  is  its  bandwidth.  When  (2.14)  is  met,  the  signal  is  considered  nar¬ 
rowband,  and  the  time  scale  compression  or  dilation  of  the  received  signal  envelope  can 
be  approximated  by  a  frequency  shift.  The  narrowband  Doppler  model  assumes,  the 
scatterer  velocity  is  small  enough,  so  that  the  entire  bandwidth  experiences  essentially 
the  same  Doppler  shift.  The  narrowband  Doppler  model  of  the  received  echo  signal  is 
then 


s{t)  «  ^SJit  -  t)  exp{j2ir(f>t) 


(2.15) 


9 


[1,  15,  22,  35,  38].  This  method  is  less  computationally  intensive  them  the  wideband 
model,  but  if  the  narrowband  criteria  are  not  met,  it  can  yield  incorrect  velocity  infor¬ 
mation  and  decreased  detection  performance  [25]. 


2.2  Echo  Processing 


The  effects  of  the  scatterer’s  range  from  the  sonar  platform  and  line-of-sight  ve¬ 
locity  on  an  echo  signal  were  discussed  in  section  2.1.  Now,  methods  of  extracting  this 
information  from  the  echo  signal  are  presented. 

A  matched  filter  is  the  optimum  maximum  likelihood  detector  for  a  slowly  fluctu¬ 
ating  point  scatterer  in  white  Gaussian  noise.  Although  it  is  optimum  only  under  these 
conditions,  it  is  a  robust  process  because  it  can  yield  useful  results  for  less  optimum  con¬ 
ditions  [3,  26,  36,  35,  39].  The  matched  filter  or  the  replica  correlator  for  the  narrowband 
and  wideband  Doppler  models  are 


nb 


2 


00 

J  r{t) f*  {t  —  f)  exp{-j2n$t)dt 
-00 


2 


^wb 

!I 

00 

J  r{t)f*  {t  -  f)]  dt 

-00 

2 


(  Narrowband  ) 


(  Wideband  ) 


(2.16) 

(2.17) 


[25,  35,  39].  The  matched  filter  is  implemented  as  the  correlation  of  the  received  signal 
with  a  bank  of  time  delayed  and  Doppler  shifted  replicas  of  the  transmitted  signal,  where 
f,  ^  and  /3  are  the  delay,  Doppler  and  dilation  factor  hypotheses  [12].  The  points  of  high 
correlation  yield  peaks  in  the  phase  plane  (delay-Doppler  plane).  This  indicates,  that 
the  received  signal  closely  matches  the  hypothesized  delay  and  Doppler  of  the  replica 
[15,  35].  furthermore,  the  shape  of  this  peak  is  determined  by  the  transmitted  waveform, 
which  ia  discussed  in  more  detail  later  in  this  section. 

The  detection  criteria  of  the  matched  filter  are 


\i\^  >  j  ^  (2.18) 

\i\^  <7  Hq,  (2.19) 


2 

where  |f|  can  be  either 
the  detection  threshold  [: 


^nb 

2 

or 

^wb 

11,  20,  33]. 

^1 

'.  and  are  detection  hypotheses  and  7  is 
1,  20,  33].  is  the  hypothesis,  that  the  sonar  signal  is  present 
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in  the  received  data,  and  is  the  hypothesis,  that  only  noise  is  present.  If  \tf  >  7, 

2 

then  h3T)Othesis  is  declared,  and  if  |^|  <  7,  then  hypothesis  is  declared. 

According  to  the  detection  hypotheses,  the  two  possible  received  signals  are 


r{t)  =  s{t)+n{t) 

r{t)  =  n(t) 


for 

for 


[15,  11,  24,  35].  The  expected  matched  filter  output  for  hypothesis  for  the  wideband 
Doppler  model  (2.17)  is 


(2.20) 


I  is{t)  +  n{t))fUp,t)dt 
-00 

where  p  =  Since  E  |ki^|  =  ^  [20],  (2.20)  becomes 


E{Kb\^\Hi}  =  El  I  {s{t)  +  n{t))f*{p,t)dt  I  {s{u)  +  niu))f*{p,u)du 


s{t)  and  n{t)  are  uncorrelated,  so  the  cross  terms  cancel,  yielding 


l~00  —00  / 

^  00  00  \ 

J  j  n{t)n*{u)f*  {p,t)  f  {p,u)dudt\  .{2.21) 

I— 00  —00  ) 


The  correlation  of  white  Gaussian  noise  with  a  power  spectral  density  of  is 


E  |n(t)n*(u)|  =  N^S  {t  —  u) 


(2.22) 


11 


[11,  35].  Therefore,  the  expected  matched  filter  output  for  an  echo  signal  backscattered 
from  a  point  scatterer  is 


1 

r 

00 

1  f 

00 

2>| 

^wb 

=  H 

1 

1 

1  s{t)f*(p,t)dt 

1 

1 

r  ^  ^  1 

1  o 

/  f*{p,t)dt 

1 

1 

-00 

J 

-00 

j 

=  E< 


00 


1-00 


1 

r 

00 

2' 

\n 

J  f*{P,t)dt 

> 

1 

1 

-OO 

d 

(2.23) 


The  noise  term  reduces  to  N^,  because  the  waveform  is  normalized  to  unity  [40],  and 
(2.23)  becomes 


E\\e. 


wb\ 


|Hl}  = 


E 


oo 


1-00 


+  N. 


(2.24) 


[22,  35],  where  E  is  the  expected  value  of  the  received  signal  energy, 
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E 

s 


(2.25) 


When  a  spread  scatterer  is  assumed  instead  of  a  point  scatterer,  (2.24)  takes  on  the  more 
general  form 


oo  00  r  °° 


-00-00 


LI-OO 


dr  d0  +  N^, 


(2.26) 


where  {p,  t)  =  .E  { (y3,  t) |^ |  is  the  scattering  function  for  a  spread  scatterer  [26]. 

Knight  et  al.  [15]  recommend  considering  the  matched  filter  process  in  the  absence 
of  noise,  because  the  properties  of  the  matched  filter  depend  upon  the  trjinsmitted  signal. 
The  output  of  the  wideband  matched  filter  normalized  to  unit  energy  in  the  absence  of 
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noise  is  known  as  the  wideband  ambiguity  function  [5,  15],  and  its  standard  form  is 

12 


00 

J  /(«)/*  [P  («  -  -r)]  du 


h-OO 


(2.27) 


(2.27)  is  the  correlation  of  the  time  dilated  and  delayed  waveform  and  the  bank  of  replicas 
[1,  25,  31,  40,  38].  Equation  (2.24)  can  be  rewritten  as  a  function  of  the  standard  form 
of  the  wideband  ambiguity  function  by  changing  the  variables  to  u  =  0{t  -  t),  0  =  ^ 
and  f  =  0{t  -  f).  Then,  (2.24)  becomes 


(2.28) 


for  a  point  scatterer.  The  narrowband  case  ^^l})  same  form  as  the 

wideband  output  except  that  the  narxowbeind  ambiguity  function 


oo 

j  (*  -  ^)  exp{-j2iT$t)dt 


-00 


(2.29) 


is  implemented  [1,  5,  15,  26,  37,  38].  When  only  the  noise  is  present  (2.19),  the  expected 
output  of  the  matched  filter  is 


{00 

0  I  n{t)f*[0{t-f)\dt 


V  1-00 


=  ^o- 


(2.30) 


as  seen  in  (2.24).  The  signal-to-interference  ratio  at  the  matched  filter  output 
can  then  be  calculated  using  (2.28)  and  (2.30), 


SIR 


wb 

<N 

^  A 

wb 

N 


(2.31) 


which  is  also  known  as  the  detection  index  (5,  15,  26,  35].  The  probability  of  detection 
for  a  fixed  probability  of  false  alarm  is 


P  =  pt+aiRo 
^  f 


(2.32) 
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[26,  35].  Therefore,  a  high  detection  index  typically  results  in  a  strong  detection  and  a 
high  level  of  confidence  in  the  estimated  scatterer  parameters. 

(2.28)  is  a  convolutional  relationship  between  the  scattering  function  and  ambigu¬ 
ity  function  in  the  delay-Doppler  plane,  which  can  be  seen  when  the  expected  matched 
filter  response  from  a  spread  scatterer  (2.26)  is  written  in  terms  of  the  ambiguity  function 
(2.27): 


E 


00  00 


//«  ^(^,T)$^^  (^0,f)drd0  +  N^ 


-00-00 


(2.33) 


[10,  23].  Therefore,  the  peak  in  the  matched  filter  response  indicates  the  delay-Doppler 
location  of  the  overlap  between  the  scattering  function  and  the  ambiguity  function  of  the 
transmitted  signal.  The  expected  matched  filter  output  is  at  a  maximum,  when  the  delay 
and  dilation  of  the  echo  signal  matches  the  replica  (^0,  f  =  ^,  ,  and  (^E 1 \H^  | ^ 
becomes 


E 


=  E  +N  , 

s  o 


(2.34) 


which  yields  a  peak  in  the  phase  plane  at  {P,t).  This  also  illustrates  the  importance 
of  the  ambiguity  function  in  the  matched  filter  response,  from  which  the  scatterer’s 
parameters  are  estimated. 

The  matched  filter  response  depends  heavily  upon  the  ambiguity  function,  as  seen 
in  (2.28).  The  estimated  scatterer  parameters  are  defined  by  the  delay-Doppler  location 
of  the  matched  filter  peak.  The  peak  covers  a  range  of  delay  and  Doppler  values,  within 
which  the  ambiguity  function  still  yields  a  significant  magnitude  relative  to  the  maximum 
response,  as  seen  in  figure  2.1  [32].  This  area  in  the  phase  plane  is  known  as  the  mainlobe, 
and  it  represents  the  uncertainty  of  the  delay  and  Doppler  estimates  for  the  scatterer. 
The  mainlobe  region  is  defined  as  the  area  within  the  half  power  points  of  ((/>,  r), 
which  is  the  delay-Doppler  cross-section  of  the  ambiguity  function  3dB  down  from  the 
peak,  and  this  is  the  resolution  area  of  the  ambiguity  function  for  a  given  waveform 
[4,  32].  Therefore,  the  resolution  of  the  delay  and  Doppler  estimates  depends  upon 
the  shape  of  the  matched  filter  peak,  and  the  shape  of  the  peak  is  determined  by  the 
ambiguity  function,  which  is  unique  to  each  sonar  signal  [5,  32,  37]. 
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Fig.  2.1.  Example  of  the  Region  of  Signal  Ambiguity  of  an  Ambiguity  Function  (dotted 
line  indicates  the  ambiguity  function  at  threshold) 


Ideally,  a  thumbtack  ambiguity  function  that  approaches  a  delta  function  yields 
more  accurate  parameter  estimations  [1,  10,  22,  35],  because  narrower  ambiguity  function 
peaks  provide  more  localized  delay  and  Doppler  estimates.  However,  ambiguity  functions 
exhibit  a  conservation  of  volume  property,  which  states  that  narrowband  signals  with 
unity  normalized  energy  have  ambiguity  functions  with  unity  volume  [5,  15,  26,  32,  37]. 
Therefore,  narrowing  the  mainlobe  results  in  higher  sidelobe  levels,  which  can  result  in 
several  sidelobes  above  the  detection  threshold,  for  some  waveforms.  The  advantages 
and  disadvantages  of  different  ambiguity  functions  are  discussed  in  section  2.3. 

2.3  General  Waveform  and  Ambiguity  Function  Characteristics 

A  wide  variety  of  waveforms  are  implemented  in  active  sonar.  For  the  purpose 
of  this  thesis,  three  types  of  sonar  waveforms  are  considered;  the  continuous  waveform 
(CW),  the  linear  frequency  modulated  waveform  (LFM)  and  the  biphase  shift  keyed 
waveform  (BPSK).  Each  individual  waveform  has  specific  attributes,  that  make  them 
advantageous  for  producing  particular  echo  information.  These  attributes  are  demon¬ 
strated  by  examining  the  ambiguity  functions  of  the  different  sonar  waveforms. 

The  discrimination  of  scatterers  that  are  separated  in  delay  and  Doppler  depends 
upon  the  resolution  of  the  sonjir  waveforms  [5,  15].  The  narrowband  Doppler  resolution 
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is 


and  the  delay  resolution  is 


Ar 

r 


(2.35) 


(2.36) 


[10,  15,  32,  35,  37].  The  accuracy  of  delay  and  Doppler  estimates  gathered  with  a 
matched  filter  is  bounded  by  the  resolution  of  a  given  sonar  signal.  Different  signals  are 
better  at  discriminating  between  delay  and  Doppler  due  to  their  length  in  time  and  their 
bandwidth  [1],  but  the  shape  of  the  peak  in  the  ambiguity  function  is  more  crucial  for 
making  unambiguous  measurements,  as  discussed  in  section  2.2. 

The  continuous  waveform  (CW)  or  sinewave  (2.1)  has  the  form 


fit)  =  a{t)  expiju^t),  (2.37) 

when  ip{t)  =  0  [15,  30].  Its  time  series  and  frequency  spectrum  are  illustrated  in  figures 
2.2a  and  b.  Sinewaves  have  a  time  bandwidth  product  of  unity  {TW  =  1),  which 


Fig.  2.2.  Example  CW  with  Center  Frequency  of  500Hz:  [a]  Time  Series  and  [b] 
Frequency  Spectrum 
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influences  their  ambiguity  functions  (figures  2.3  and  2.4)  [10,  32].  The  ambiguity  function 


Fig.  2.3.  Long  CW  Ambiguity  Function  Fig.  2.4.  Short  CW  Ambiguity  Function 


Fig.  2.5.  Anibiguity  Function  Cross-Sections:  [a]  Long  CW,  [b]  Short  CW 


peak  for  a  sinewave  is  spread  in  delay  for  the  time  duration  of  the  signal  (T)  and  spread  in 
Doppler  according  to  the  bandwidth  of  the  signal  (^W  =  This  results  in  a  CW  with 
a  large  T  correlating  over  a  large  range  of  delays  but  only  a  small  range  of  frequencies, 
as  seen  in  figure  2.5a  [4,  32],  If  the  CW  has  a  relatively  short  T,  then  it  correlates  over 
only  a  small  range  of  delays  but  over  a  large  range  of  frequencies,  as  seen  in  figure  2.5b. 
The  CW  can  have  good  Doppler  resolution  with  a  longer  CW  or  good  range  resolution 
with  a  shorter  CW,  but  not  both  simultaneously  [15]. 
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and  typically  TW  >  1  [15,  32].  Since  the  frequency  of  the  LFM  increases  linearly  with 
time,  the  LFM  correlates  over  short  time  intervals  at  different  frequency  intervals,  which 
result  in  a  diagonal  ambiguity  function  peak  across  the  phase  plane  with  an  overall 
mainlobe  frequency  range  of  W  and  a  delay  range  of  T,  as  seen  in  figures  2.7  and 
2.8  [4,  5,  23,  32].  The  LFM  is  not  a  good  joint  estimator  of  both  delay  and  Doppler 
parameters,  because  it  is  highly  skewed  in  the  delay-Doppler  plane.  The  minimum 
variance  of  the  ambiguity  function  estimation  in  delay  and  Doppler  depends  upon 

$^^(/3,-r)  =  $^^(/3,T),  (2.39) 

which  is  satisfied  by  even  symmetric  waveform  functions  [21].  When  (2.39)  is  satisfied, 
the  waveform  is  said  to  be  uncoupled,  and  the  variance  of  the  delay  and  Doppler  estima¬ 
tions  are  minimaJ.  Otherwise,  the  waveform  is  coupled,  and  this  coupling  is  represented 
by  the  slope  in  the  ambiguity  function  [7].  Therefore,  the  LFM  is  a  highly  coupled  wave¬ 
form,  so  the  joint  delay  and  Doppler  estimates  for  a  LFM  can  not  be  localized  better 
than  T  and  W,  if  both  of  the  parameters  are  unknown  [5,  35].  However,  if  one  of  the 
parameters  is  already  known  (0  or  t),  the  other  can  be  estimated  with  a  much  higher 
resolution  (^  or  ^)  [32]. 

The  BPSK  is  a  combination  of  CW  chips  at  the  same  frequency,  that  are  phase 
shifted  by  ±180°  according  to  a  pseudo- random  code  in  order  to  create  a  broadband 
waveform  that  appears  noise- like,  as  seen  in  figures  2.9a  and  b  [22].  The  3dB  bandwidth 
of  the  BPSK  is  W  «  where  N  is  the  length  of  the  pseudo-random  code  sequence 
(or  number  of  chips).  An  M  stage  shift  register  generates  the  code  sequence  of  length 
N  =  2^  —  1,  that  shifts  between  ±180°  and  repeats  for  eew;h  chip  [16,  22].  The  frequency 
spectrum  of  a  BPSK  waveform  appeairs  similar  to  a  sine  function  (sinc(x)  =  sin(i)/ar), 
where  the  mainlobe  has  the  bandwidth  of  W,  and  the  first  set  of  sidelobes  are  approxi¬ 
mately  13dB  down,  as  seen  in  figure  2.9b.  The  BPSK  has  very  good  joint  frequency  and 
delay  resolution  (^,  because  it’s  ambiguity  function  has  an  exceptionally  localized 
mciinlobe  (figures  2.10  and  2.11).  Furthermore,  it  is  uncoupled  in  delay  and  Doppler, 
when  the  complex  spectrum  of  the  BPSK  is  even  and  satisfies  (2.39).  The  phase  plane 
area  of  the  mainlobe  of  the  ambiguity  function  for  the  BPSK  is  which  can  be  very 
high  resolution  [4,  32].  This  higher  resolution  can  make  it  more  computationally  inten¬ 
sive  to  search  the  phase  plane  for  detection  peaks,  and  it  results  in  a  very  large  sidelobe 
region  for  the  BPSK. 
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Fig.  2.9.  Example  BPSK  consisting  of  30  Chips  with  Center  Frequency  of  IkHz  and  a  Bandwidth 
of  lOOHz  :  [a]  Time  Series  and  [b]  Frequency  Spectrum 


Fig.  2.10.  BPSK  Ambiguity  Function  Fig.  2.11.  BPSK  Ambiguity  Function 

Cross-Section  (solid  line  is  the  mainlobe  re¬ 
gion,  dashed  line  is  the  sidelobe  region) 
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The  detection  of  a  scatterer  in  reverberation  and  the  detection  of  multiple  scat- 
terers  closely  spaced  in  delay  and  Doppler  (or  a  spread  scatterer)  are  two  particular  cir¬ 
cumstances,  that  illustrate  some  of  the  advantages  and  disadvantages  of  different  sonar 
signals.  Reverberation  is  backscattered  energy  from  discontinuities,  boundaries  and  the 
volume  in  the  ocean  [3,  35].  These  discontinuities  range  from  water-born  particulates  to 
biological  sources  and  from  surface  interactions  to  bottom  backseat tering.  Reverberation 
is  represented  by  the  scattering  function 

R^{t,4>)  =  e{\s^{t,<I>)\^]  (2.40) 

[23].  5^(r,  0)  is  the  reverberation  spreading  function,  which  is  similar  to  (2.25),  ex¬ 
cept  that  it  accounts  for  the  spread  scattering  from  all  of  the  ocean  inhomogeneities 
encountered  instead  of  just  a  point  scatterer.  The  expected  matched  filter  response  for 
reverberation  is 

00  00 

-00  -00 
00  00 

=  I  ^{|5^(r,0)p}$^^(0,f)d<^dT,  (2.41) 

-OO  — OO 

which  is  the  same  form  as  the  scattering  from  a  spread  scatterer  (2.33)  [23,  35].  Rever¬ 
beration  is  the  dominant  form  of  interference  at  short  ranges  from  the  sonar  platform, 
but  its  magnitude  decays  with  range,  like  all  backscattered  signals  [36].  Once  a  range 
has  been  reached,  where  reverberation  has  decayed  below  the  ambient  noise  level,  noise 
becomes  the  dominant  form  of  interference  [23].  Reverberation  occurs  at  low  Doppler, 
because  the  discontinuities  in  the  ocean  aire  stationary  or  moving  at  very  low  velocities. 
Therefore,  reverberation  backscatter  is  confined  to  a  band  of  width  B  around  the  trans- 
mit  frequency  over  the  entire  receive  interval  in  the  phase  plane,  as  seen  in  figures  2.12a 
and  b  [4,  23].  This  Doppler  bandwidth  is  determined  by  the  motion  of  the  medium, 
and  it  can  be  modeled  for  the  purpose  of  performance  prediction  as  uniform  over  the 
entire  bandwidth  [23].  Since  reverberation  is  often  the  main  source  of  interference  in  ac¬ 
tive  sonar  [36],  its  effects  need  to  be  considered,  when  evaluating  the  utility  of  different 
sonar  signals. 
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Fig.  2.12.  Matched  Filter  Responses  of  Different  Signals  in  Reverberation:  [a]  low 
Doppler  scatterer,  [b]  high  Doppler  scatterer 


High  levels  of  reverberation  can  become  prominent  in  matched  filter  responses, 
as  seen  in  figures  2.12a  and  b.  The  degradation  of  the  matched  filter  response  due  to 
reverberation  is  determined  by  the  amount  of  overlap  between  the  ambiguity  function 
of  the  signed  and  the  reverberation  scattering  function  [4,  18,  35].  Therefore,  a  sonar 
with  minimal  overlap  between  (j>,  and  (r,  0)  in  a  given  scenario  is  chosen  to 
majcimize  detection  [15].  The  ambiguity  function  mainlobe  is  embedded  in  reverberation 
when  the  scatterer  has  a  low  velocity.  The  ambiguity  function  mainlobes  of  short  CW 
or  LFM  signals  are  spread  in  Doppler  to  a  length  of  W,  which  may  not  be  completely 
embedded  in  reverberation,  it  <  W  (figure  2.12a).  The  ambiguity  functions  of  the 
short  duration  CW  and  LFM  have  less  overlap  with  the  reverberation  scattering  function 
than  the  long  CW  and  BPSK,  due  to  their  Doppler  spread  mainlobes.  Therefore,  the 
effects  of  reverberation  are  minimized,  and  there  is  a  greater  chance  that  the  scatterer 
parameters  can  be  estimated  from  them  [4,  15,  18,  35].  However,  long  CWs  and  BPSKs 
may  provide  better  detection  in  high  Doppler  sceneirios,  because  their  mainlobes  are 
less  spread  in  Doppler  (figure  2.12b).  The  long  CW  may  have  a  sufficiently  narrow 
bandwidth,  such  that  its  ambiguity  function  has  negligible  overlap  with  the  reverberation 
scattering  function  at  high  Doppler  [4,  18].  The  mainlobe  of  the  BPSK  may  also  have 
minimal  overlap  with  the  reverberation  scattering  function,  but  its  Ijirge  sidelobe  region 
will  most  likely  overlap  [15,  18,  35].  This  makes  the  BPSK  a  less  effective  waveform  for 
the  detection  of  a  high  velocity  scatterer  than  a  long  duration  CW.  However,  it  may  still 
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be  more  effective  then  the  LFM  or  the  short  CW,  because  its  mainlobe  can  have  less 
overlap  with  the  reverberation  scattering  function. 

The  discrimination  of  multiple  scatterers  that  are  closely  spaced  in  delay  and 
Doppler  is  often  of  interest  in  active  sonar.  Signals  with  fine  time  resolution  are  needed 
to  discriminate  between  multiple  scatterers  closely  spaced  in  delay,  and  signals  with  high 
frequency  resolution  discriminate  between  multiple  scatterers  closely  spaced  in  Doppler, 
cis  seen  in  (2.35)  and  (2.36).  Clearly,  a  signal  with  good  resolution  in  time  and  frequency 
is  required  to  successfully  discriminate  individual  scatterers  closely  spaced  in  delay  and 
Doppler.  For  this  reason,  BPSKs  are  more  appropriate  than  CW  or  LFM  signals  for  the 
discrimination  of  multiple,  closely  spaced  scatterers. 

2.4  Motivation  behind  Sub-band  Modulation 

Careful  waveform  selection  is  required  in  order  to  accurately  estimate  the  scat- 
terer  parameters  under  particular  conditions.  Therefore,  different  waveforms  are  needed 
to  perform  detection  and  estimation.  It  can  be  beneficial  to  transmit  numerous  active 
sonar  signals  at  once,  in  order  to  gather  more  information  about  a  scatterer.  Numerous 
variations  of  signal  combinations  could  be  used  to  gather  different  types  of  echo  informa¬ 
tion.  An  obvious  choice  is  to  transmit  signals  with  good  frequency  resolution  and  signals 
with  good  time  resolution  together,  in  order  to  gather  high  resolution  delay  and  Doppler 
information.  Another  alternative  is  to  transmit  a  high  resolution  and  low  resolution  sig¬ 
nal  together,  in  order  to  decrease  detection  complexity.  A  low  resolution  signal  will  result 
in  a  spread  peak  in  the  matched  filter  response,  that  does  not  localize  well.  The  rough 
localization  of  delay  and  Doppler  parameters  of  the  scatterer  can  be  made  from  the  low 
resolution  response,  and  applied  to  the  match  filtering  of  the  high  resolution  signal.  The 
computational  requirements  to  secirch  a  large  delay-Doppler  region  with  a  bank  of  high 
resolution  matched  filter  replicas  are  considerably  reduced  by  first  isolating  the  scatter 
with  a  low  resolution  matched  filter.  Knight  [15]  noted  that,  a  combination  of  different 
sonar  signals  are  usually  transmitted  to  detect  tcirgets  of  various  velocities.  Therefore, 
a  signal  with  good  reverberation  performance  but  low  Doppler  resolution  and  a  signal 
with  good  Doppler  resolution  but  poor  reverberation  performance  is  another  possible 
signal  combination.  This  signal  combination  allows  for  a  detection  to  still  be  made  in 
a  high  reverberation/low  Doppler  scenario,  as  well  as  achieve  good  Doppler  information 
with  a  higher  scatterer  velocity.  The  desire  to  transmit  multiple  sonar  signals  together 
is  one  of  the  motivating  factors  behind  sub-band  modulation. 
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Transmission  time  is  limited  in  monostatic  sonar.  There  must  be  separate  trans¬ 
mit  and  receive  intervals,  because  the  trcinsmit  and  receive  transducers  are  colocated 
(and  are  most  likely  the  .same  transducers).  Scatterers  within  the  blind  range 


will  not  be  detected,  since  signals  can  not  be  received  while  transmitting.  The  transmit 
interval  is  limited  to  minimize  the  blind  range  [30,  15].  However,  it  is  beneficial  to  utilize 
the  longest  transmit  time  possible,  because  frequency  resolution  is  proportional  to  the 
inverse  of  time  (2.35).  This  presents  a  problem,  when  transmitting  consecutive  sonar 
signals  during  the  same  transmit  interval,  because  the  transmit  time  of  each  signal  is 
limited  even  further.  Maximizing  the  transmission  time  of  each  sonar  signal  is  another 
motivation  behind  sub-band  modulation. 

One  of  the  purposes  of  transmitting  multiple  sonar  signals  is  to  receive  more  echo 
information  describing  a  specific  event;  for  example,  a  scatterer’s  location  and  velocity 
at  a  given  moment  in  time.  The  gathering  of  echo  information  with  time  separated 
sonar  signals  can  differ  between  each  interrogation  due  to  the  nonstationarity  of  echo 
parameters.  The  time  separation  of  the  soneir  signals  and  the  potential  for  a  moving 
scatterer  can  result  in  echo  information  describing  multiple  events,  as  opposed  to  a 
specific  event.  Clearly,  the  most  beneficial  method  of  achieving  this  goal  is  through 
simultaneous  signed  transmission,  which  also  provides  the  maximum  transmission  time 
for  each  sonar  signal. 

Transmit  energy  is  another  concern  in  active  sonar.  The  transmit  energy  of  the 
waveform  is  defined  as 

T  T  T 

=  Jp{t)dt  =  |[f^{t)fdt  =  (2.43) 

0  0  0 

where  P{t)  is  the  instantaneous  transmitted  power  [11,  30].  However,  is  energy 
only  when  the  waveform  is  normalized  (2.3).  The  ^  relates  to  the  peak  of  the  real 
transmitted  waveform,  and  is  dropped  when  the  the  real  waveform  is  normalized  instead 
of  the  analytic  waveform.  A  maximum  operational  level  is  set  for  each  peak  power 
limited,  sonar  projector.  If  this  maximum  level  is  surpassed,  cavitation  may  occur. 
Cavitation  is  the  formation  and  collapse  of  bubbles  on  the  face  of  the  sonar  projector, 
which  results  in  the  deterioration  of  projector  performance  and  an  impedance  mismatch 
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between  the  projector  and  the  medium  around  it  [3,  36].  Detection  is  directly  related  to 
the  received  signal-to-interference  ratio  SIR^,  which  in  turn  is  related  to  the  attenuated 
transmit  level  [4,  36].  The  transmit  energy  needs  to  be  maximized  to  produce  a  higher 
probability  that  the  sonar  signal  will  be  detected  in  noise  limited  conditions.  The  overall 
transmit  energy  can  be  related  to  this  maximum  transmit  level  through  energy  efficiency 
.  Energy  eflficiency  is  the  transmit  energy  of  a  real  waveform  compared  to 
the  maximum  transmit  energy 
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(2.44) 


where  /  is  the  peak  amplitude  of  the  real  transmit  waveform.  Therefore,  energy 

7TICLX 

efficient  is 

^ef,  = 

peak  max 


where  is  the  peak  transmitted  power  [30].  Energy  efficiency  is  maximized  in  the 

case  of  a  single  real  sinusoidal  signal  normalized  to  where 


E^  e 

E  ..  =  -=r-^  =  -2  =  1.  (2.46) 

effEe  ^  ' 

max  o 


When  transmit  energy  is  maximized,  the  probability  of  detection  for  a  fixed  probability 
of  false  alarm  is  also  maximized  under  noise  limited  conditions,  as  seen  in  (2.32).  How¬ 
ever,  the  total  transmit  energy  is  not  always  maximized,  when  multiple  sonar  waveforms 
axe  simultaneously  in  a  peak  power  limited  system.  This  can  be  caused  by  amplitude 
modulation,  and  is  discussed  in  section  3.1.  Therefore,  another  motivational  goal  behind 
sub-band  modulation  is  to  avoid  amplitude  modulation.  In  summary,  this  motivation 
conveys  the  need  to  investigate  a  process,  that  permits  the  time  simultaneous  transmis¬ 
sion  of  multiple  soneir  signals,  while  avoiding  amplitude  modulation. 
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Chapter  3 

Modulation  Theory 

Modulating  a  carrier  waveform  with  several  active  sonar  waveforms  is  one  way  to 
transmit  multiple  waveforms  simultaneously.  This  has  been  done  in  radio  broadcasting 
for  many  years.  Both  amplitude  and  angle  modulation  theory  are  discussed  and  the 
sub-band  modulation  concept  is  defined  and  analyzed.  This  is  followed  by  an  in-depth 
analysis  of  sub-band  modulation.  Finally,  the  demodulation  of  sub-band  modulated 
signals  is  thoroughly  examined. 

3.1  Overview  of  Amplitude  Modulation 

Linear  combination  is  one  way  to  simultaneously  transmit  multiple  sonar  wave¬ 
forms  without  sacrificing  the  transmit  time  of  each  individual  waveforms.  This  is  done 
by  selecting  waveforms,  that  are  sufliciently  spectrally  separated  to  permit  their  recovery 
by  filtering.  The  linear  combination  of  multiple,  analytic  sonar  waveforms  is 

N  _  .  . 

fi*)  =  <e(0,T) 

m=l  \  m  / 

N  .  . 

^  S  0  <e(0,r),  (3.1) 

m=l  ^  ^ 

where  /^(t)  are  the  complex  envelopes  of  the  sonar  waveforms,  and  N  is  the  total  number 
of  sonar  waveforms.  The  transmit  waveform  energy  is  normalized  by 


since  the  real  waveform  will  be  trcinsmitted,  and  the  transmit  energy  is  limited  to 
[30].  depends  on  the  modulating  signals  and  modulation  type,  so  it  will  be  redefined 
for  each  form  of  amplitude  modulation  discussed  in  this  section.  (3.1)  can  be  rewritten 
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as 


r  1 

\ 

m=l  ^ 

/ 

CJ  — 

0 

L  m 

0 

(3.3) 


by  separating  the  carrier  waveform.  (3.3)  reveals,  that  the  linear  combination  of  sonar 
waveforms  results  in  a  carrier  waveform  with  a  modulating  amplitude.  Although  one 
of  the  motivations  behind  sub-band  modulation  is  to  avoid  amplitude  modulation,  it  is 
presented  here  to  illustrate  why  it  is  undesirable. 

If  the  modulation  is  a  real  function  instead  of  analytic,  then  (3.3)  becomes  am¬ 
plitude  modulation  (AM).  Amplitude  modulation  is  the  linear  combination  of  multiple 
waveforms  as  the  amplitude  component  of  a  carrier,  and  has  the  form 


>/2 

/(<)  =  ye;A^exp(;a;o<) 

771=1 

(3.4) 


where  w  =  ui  -  J  >  w.  and  N  is  the  total  number  of  sonar  waveforms  in  the 

771  0  U 

77Z  771  -  _ 

modulating  spectrum,  which  is  an  even  number  for  AM  suppressed  carrier  [19,  30,  41]. 
Since  the  maximum  transmit  energy  is  limited  to  (section  2.3),  from  (3.2)  becomes 


An  AM  waveform  is  typically  designed  using  (3.4)  and  not  (3.1),  so  the  u  are  the  center 
frequencies  of  the  sonar  waveforms  and  The  carrier  waveform  is  suppressed  in 

this  case  (commonly  known  as  AM  suppressed  carrier),  and  the  time  series  and  frequency 
spectrum  of  the  amplitude  modulation  of  a  suppressed  carrier  by  two  sine  waves  is  shown 
in  figures  3.1a  and  b^. 

^The  non-symmetries  in  any  of  the  frequency  spectrum  plots  in  this  thesis  are  due  to  finite 
sampling. 
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Fig.  3.1.  Example  AM  Suppressed  Carrier  of  Two  Sinewaves  (Center  Frequencies  of 
500Hz,  800Hz)  with  a  Carrier  Frequency  of  2kHz:  [a]  Time  Series  and  [b]  Frequency 
Spectrum 


Assuming  that  f^{t)  are  real  and  sinusoidal,  the  total  transmit  energy  (2.43)  of 


a  real  AM  waveform  with  a  suppressed  carrier  (3.4)  is 


e  A 


2  T 


'N/2 


12 


m=l 


0 

The  peak  amplitude  of  the  real  transmit  waveform  is 

/  =  maximum  f(t)  =  ./e~ A 

■’max  ^  '  V  o  I 


dt. 


N/2 

E 

771=  1 


producing  a  maximum  transmit  energy  (2.44)  of 


e  J 


E 


TTtCLX 


'n/2 

E 

m=l 


dt  =  e 


(3.6) 


(3.7) 


(3.8) 
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(3.6),  (3.5)  and  (3.8)  are  combined  to  yield  the  energy  efficiency  of  an  AM  suppressed 
carrier  waveform, 


^eff 


E 


0,  {t)cos(u  t  +  tp  (t)) 

|2 

dt 

1  ^ 

if 

0 

2 

dt 

T 

/ 

0 


1 2 

a  (t)cosfc*;  t  +  ip  (t)')  dt 


T 

f 

0 


dt 


(3.9) 


Since  cosine  functions  are  limited  to  a  range  of  (-1,1),  the  numerator  of  (3.9)  is  always 
less  than  or  equcil  to  the  denominator.  Therefore,  E^^^  <  1  for  AM  suppressed  carrier 
waveforms.  This  illustrates  that  the  transmit  energy  of  AM  suppressed  carrier  waveforms 
depend  on  the  modulating  waveforms,  and  their  overall  transmit  energy  is  not  maximized. 

Another  form  of  amplitude  modulation  occurs,  when  the  carrier  is  not  suppressed. 
The  AM  waveform  with  a  carrier  is 


f(t)  =  y^A^exp(ju;ot) 


(iV-l)/2 

1+  E  < 

m=l 


m 


(3.10) 


where  the  carrier  appears  in  the  spectrum  (see  figure  3.2b),  and  N  is  odd  [2,  9, 13, 19,  41]. 

A  is  redefined  for  an  AM  waveform  with  carrier  as 
c 


A  = 
c 


(3.11) 


The  totzd  transmit  energy  of  a  real  AM  waveform  with  a  carrier  (3.10)  is 


(Ar-i)/2 

0 

771  =  1 

(3.12) 
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Fig.  3.2.  Example  AM  with  Carrier  of  Two  Sinewaves  (Center  Frequencies  of  500Hz, 
800Hz)  with  a  Carrier  Frequency  of  2kHz:  [a]  Time  Series  and  [b]  Frequency  Spectrum 


By  taking  into  accoimt  the  peak  amplitude  of  the  real  transmit  waveform  of  the  AM 
waveform  with  carrier 


/  =  ,/e“A 

‘‘max  V  o  c 


iN-l)/2 

^  +  I] 

m=l 


>  /(*). 


(3.13) 


which  yields  a  maximum  transmit  energy  of 


E 


max 


e  T 
0 


(iV-l)/2 

1+  E 


m=l 


dt  =  e  . 
o’ 


(3.14) 


the  energy  efficiency  of  an  AM  waveform  with  a  carrier  can  be  calculated. 


E 


dt 


eff 


12 


dt 


T 

f 

0 


Tr 

f 

0 


l  +  ^(iV-l)/2  (i) 

m=l  ^ 


dt 


(3.15) 
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where  ^  1  under  the  same  conditions  that  applied  to  the  AM  suppressed  carrier 
case.  The  influence  of  the  modulating  waveforms  is  still  present,  and  the  overall  trans¬ 
mit  energy  is  not  maximized.  Therefore,  the  varying  amplitude  of  an  AM  waveform  is 
disadvantageous  in  a  peak  power  limited  system. 


3.2  Overview  of  Angle  Modulation  and  Sub-band  Modulation 


Multiple  sonar  waveforms  can  also  be  transmitted  simultaneously  using  angle 
modulation.  Wherein,  the  amplitude  of  the  carrier  waveform  is  held  constant,  and 
multiple  waveforms  are  combined  as  the  phase  component  of  the  carrier  [2,  9,  41].  The 
use  of  angle  modulation  in  active  sonar  was  first  suggested  by  L.  H.  Sibul  [30].  Angle 
modulation  may  be  either  phase  modulation  (PM)  or  frequency  modulation  (FM).  The 
phase  of  the  cjirrier  waveform  is  varied  linearly  with  the  message  waveforms  with  PM, 
where 

(3.16) 


N 

a;  t  +  k  y' 

0  p  ^ 

m=. 


is  the  deviation  constant  of  the  phase  modulation  and  will  be  further  explained  later 

on  in  this  section  [2,  27,  9,  17,  28,  34,  41].  is  defined  as  A^  =  1/V^,  because  the 

amplitude  of  an  angle  modulated  waveform  remains  constant  [9],  and  the  the  maximum 

transmit  energy  is  limited  to  e  .  Middleton  [17]  states,  that  the  conditions  u  k 

o  11^  m  p 

Wq  and  must  be  imposed,  if  (3.16)  is  to  be  characterized  as  modulation  and  not 

waveform  mixing.  An  example  of  the  phase  modulation  of  a  carrier  with  two  sinewaves 
is  shown  in  figures  3.3a  and  b. 

FM  is  the  linear  variation  of  the  time  derivative  of  the  phase  of  the  carrier,  as 
shown  by 


^{/(*)}  = 


\ 

du 

J 


(3.17) 


where  A^  remains  defined  as  k ^  is  the  deviation  constant  of  the  frequency 

modulation  [9,  28,  17,  19,  30,  41],  and  the  modulation  conditions  for  FM  are  k ^  < 
Wp  cind  0/^  Wq  [17].  An  example  of  two  sinewaves  frequency  modulating  a  carrier  is 
shown  in  figures  3.4a  and  b. 
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Fig.  3.3.  Example  PM  of  Two  Sinewaves  (Center  Frequencies  of  500Hz,  800Hz)  with 
a  Carrier  Frequency  of  2kHz  and  a  Modulation  Index  of  1.5:  [a]  Time  Series  and  [b] 
Frequency  Spectrum 


The  transmit  energy  (2.43)  for  real  waveforms  with  both  real  phase  and  frequency 
modulation  is 


1 

0 


(3.18) 


Since  the  amplitude  of  the  carrier  waveform  is  not  changed  by  the  frequency  and  phase 
modulating  waveforms,  the  transmit  energy  becomes  [9, 41].  Certain  extenuating 

circumstances  do  exist,  which  keep  from  being  maximized.  These  circumstances 
become  more  obvious,  when  (3.18)  is  expressed  as  a  series  of  Bessel  functions.  This 
will  be  discussed  thoroughly  in  section  3.2.1.  However,  they  can  be  typically  avoided  by 
limiting  the  size  of  so  the  power  of  an  augle  modulated  waveform  is  independent  of  the 
message  waveforms  [41].  Angle  modulation  transfers  the  signal  energy  from  the  carrier 
to  the  sidebands,  so  the  average  energy  of  the  modulating  waveform  is  the  same  as  the 
unmodulated  carrier  waveform  [9,  17].  Therefore,  it  is  no  surprise,  that  energy  efficiency 
for  angle  inodulation  is  ~  transmit  energy  is  maximized. 

The  difference  between  the  trcinsmit  energy  in  au  amplitude  modulated  waveform 
and  the  transmit  energy  in  an  angle  modulated  waveform  can  be  quantified  by  evalu¬ 
ating  the  contribution  of  each  individual  waveform  to  the  signal  energy.  AM  splits  the 
amplitude  of  the  ceurrier  waveform  amongst  the  waveforms  modulating  it.  When  equal 
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[a]  time  (s)  [b]  frequency  (Hz) 


Fig.  3.4.  Example  FM  of  Two  Sinewaves  (Center  Frequencies  of  500Hz,  800Hz)  with 
a  Carrier  Frequency  of  2kHz  and  a  Modulation  Index  of  0.15:  [a]  Time  Series  and  [b] 
Frequency  Spectrum 


amplitude  signals  are  used,  the  amplitude  for  e2w;h  waveform  is  Therefore,  the  signal 
energy  for  each  of  the  individual  waveforms  is  reduced  by  a  factor  of  This  differs 
from  angle  modulation,  where  the  signal  energy  instead  of  the  amplitude  is  split  between 
the  waveforms  modulating  the  carrier.  Assuming  that  all  of  the  signal  energy  is  split 
between  the  waveforms  modulating  the  carrier  and  none  is  lost,  then  energy  in  each 
waveform  becomes  -^E^,  instead  of  '^E^.  Therefore,  the  energy  splitting  of  angle  mod¬ 
ulation  provides  an  increase  in  individual  waveform  energy  over  the  amplitude  splitting 
of  amplitude  modulation. 

Angle  modulation  allows  for  the  simultaneous  transmission  of  multiple  sonar  wave¬ 
forms  over  a  maximized  time  interval  while  avoiding  amplitude  modulation.  However, 
this  research  focuses  on  phase  modulation,  rather  than  frequency  modulation,  for  several 
reasons.  Most  sonar  waveforms  are  combinations  of  sinewaves,  in  one  form  or  another. 
Whether  they  are  simple  CWs  or  a  combination  of  CW  pulses  used  to  represent  a  LFM, 
sinusoidal  waveforms  make  up  the  majority  of  the  waveforms  used  in  this  study.  The 
integration  of  a  sinewave  is  another  sinewave  phase  shifted,  so  ip  the  context  of  this 
research,  frequency  modulation  and  phase  modulation  are  essentially  the  same  process. 
However,  frequency  modulation  requires  more  computational  processing.  The  processing 
is  simplified  through  the  use  of  PM,  because  it  requires  a  less  elaborate  implementation 
in  Matlab.  Therefore,  the  definition  of  sub-band  modulation  for  the  purpose  of  this 
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thesis  is  the  phase  modulation  of  a  carrier  by  multiple  waveforms,  as  applied  to  active 
sonar. 


3.2.1  Characteristics  of  Sub-band  Modulation 

The  two  major  parameters  of  phase  modulation  as  described  in  (3.16)  are  the 
carrier  frequency  and  the  deviation  constant  k^.  The  carrier  frequency  is  the  center 
frequency  of  the  sub-band  modulated  waveform,  and  the  modulating  waveforms  appear 
in  the  frequency  spectrum  as  the  lower  and  upper  sidebands  around  the  carrier  frequency, 
as  seen  in  figure  3.5  where  —  ^  —  2kHz. 

The  deviation  constant  is  the  maximum  value  of  the  phase  deviation  of  the  mod¬ 
ulated  carrier  waveform.  When  the  message  waveforms  are  sinusoidal,  it  is  known  as 
the  modulation  index  [41],  and  in  this  thesis  the  terms  will  be  used  interchangeably. 
The  deviation  constant  controls  the  bandwidth  of  the  modulation  and  the  magnitude  of 
the  sidebands,  and  therefore  controls  the  prominence  of  the  sub-band  waveforms  in  the 
modulating  spectrum.  The  effects  of  are  apparent  when  (3.16)  is  expressed  as  the 
sum  of  Bessel  functions.  First,  (3.16)  is  rewritten  as 

I  ^ 

wii)} = n 

I  m=l 


Assuming  that  the  sonar  waveforms  f^{i)  sinusoidal,  the  Bessel  function  expansion 
is 

00 

exp  [^jfe^sin  (wt)]  =  -/^(y  exp(jma;t),  (3.19) 

m=—oo 


and  (3.16)  becomes 

I  (3.20) 


00  00 

=  E  •••  E  4 

771^  =— 00  m^  =  -00 


{k  )---J.  {k  )cos([a;„- 
\  P  P  U  0 
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nm)  = 


N 


exp 
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m=l 


00 


m=“00 


4  (V“‘’('’'™"ra‘) 

m  ^ 


[2,  9,  17,  19,  28,  41].  The  sub-band  modulating  waveforms  lie  at  the  carrier  frequency 
plus  and  minus  integer  multiples  of  their  center  frequencies  with  magnitudes  determined 
by  the  value  of  the  Bessel  function  at  the  deviation  constant.  The  form  of  the  sub-band 
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modulation  is  shown  by  the  following  example,  where  two  sinusoidal  message  waveforms 
are  sub-band  modulating  a  carrier.  (3.20)  is  rewritten  as 


OO  00 


a;o  +  mia;i  + 77120^2 


0’ 


m^  =  -oom2  =  ~oo 


(3.21) 


and  is  displayed  in  figure  3.5. 


Fig.  3.5.  Sub-band  Modulation  of  Two  Sinewaves  (r/^^  =  450Hz,  ri^  =  500Hz)  with  77^  = 
2kHz  and  =  0.7  (Frequency  Spectrum) 
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Ziemer  and  Tranter  [41]  point  out,  that  angle  modulation  results  in  an  infinite 

number  of  sidebands,  which  yields  an  essentially  infinite  bandwidth.  However,  a  fixed 

yields  lim  J  .(A:  )  =  0,  so  only  the  sidebands  containing  significant  power  are  applica- 
m->oo  m  p 

ble.  K  represents  the  significant  number  of  sidebands  on  each  side  of  the  carrier,  where 
J^k  )  76  0  for  m  G  (-«,«).  Initial  investigations  into  phase  modulation  with  Matlab 
revealed,  that  even  k  sidebands  ~40dB  down  from  the  carrier  are  significant  and  aid  in 
the  proper  demodulation  of  the  sonar  signals.  «  can  be  related  to  the  bandwidth  of  the 
sub-band  modulation  by 


B  =  2ku} 


max' 


(3.22) 


where  ui  is  the  highest  frequency  of  the  modulating  waveforms,  and  B  is  the  band- 

TfXQtX 

width  of  the  modulation  [41].  Carson’s  Rule, 


B 


-2(*  +1) 


UJ 


max 


(3.23) 


which  describes  the  bandwidth  of  angle  modulation  as  a  function  of  the  deviation  con¬ 
stant  [4,  41],  relates  k  and  k^  by 


K^k^  +  1,  (3.24) 

illustrating  that  the  bandwidth  of  the  sub-band  modulation  broadens  with  increasing 

values  of  k  . 

P 

This  can  be  problematic  in  sub-band  modulation,  because  the  transmit  energy  can 
be  effected  by  large  increases  in  modulation  index  and  bandwidth.  The  overall  transmit 
energy  of  a  sub-band  modulating  waveform  (3.18)  can  be  restated  using  (3.20)  to  better 
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illustrate  this  effect; 


(3.25) 


The  second  term  of  f  (3.25)  is  negligible  for  small  values  of  f  [29])  ,  and 

(ifc  )  w  1  [9,  41].  The  second  term  becomes  significant  as  k  becomes  larger,  which 
0  P  P 

also  results  in  a  decrease  in  the  magnitude  of  the  first  term  [17,  29].  Schwartz  [29]  defines 
a  large  modulation  index  as  k^>\.  t  can  still  approximate  unity  in  this  case.  However, 
when  the  modulation  index  becomes  too  large,  the  modulating  bandwidth  becomes  so 
wide,  that  some  of  the  significant  terms  of  (3.25)  are  located  at  frequencies  much  higher 
than  the  carrier  Wq.  This  can  result  in  f  56  1,  which  equates  to  amplitude  modulation. 
Computational  complexity  and  sampling  frequency  limitations  may  require  some  of  these 
higher  frequency  terms  to  be  neglected,  which  can  also  result  in  amplitude  modulation 
[13]. 

Therefore,  the  modulating  waveforms  begin  to  affect  the  overall  transmit  energy 

as  k  increases,  as  shown  by  the  following  example.  The  sub-band  modulation  of  a  carrier 
P 

with  two  sinewaves  (3.21)  is  illustrated,  and  the  time  series  and  spectra  of  the  sub-band 

waveform  with  increasing  k^  are  shown  in  figures  3.6a- f.  Small  values  of  k^  {k^  = 

1.2)  result  in  a  narrow  bandwidth  of  essentially  only  the  first  sideband  {B  ss  2.2kHz), 

as  seen  in  figures  3.6a,  b.  Here,  the  spectrum  of  the  sub-band  modulation  appears 

somewhat  similar  to  an  AM  spectrum,  because  most  of  the  waveform  information  is 

contained  within  the  major  sidebands.  However,  the  sidebands  quickly  broaden  and 

begin  to  resemble  overlapping  Bessel  functions  as  k^  increases  (figures  3.6c,  d).  The 

waveform  information  is  no  longer  contained  within  the  first  sidebands,  when  larger 

values  of  k  are  used  (k  =  2,  B  «  3kHz).  Therefore,  sub-band  modulation  can  be 
P  P 

used  to  create  a  waveform  with  a  spread  spectrum.  Eventually,  the  bandwidth  of  the 


Magnitude  Magnitude  Magnitude 
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sub-band  modulation  will  fill  the  spectrum  as  k  increases  further  (k  =  5.5,  B  » 

P  P 

6.5kHz),  as  seen  in  figures  3.6e  and  f.  This  can  even  result  in  a  modulating  amplitude 
in  some  extreme  cases,  as  seen  in  figure  3.6e.  When  k^  is  large,  an  angle  modulating 
waveform  can  be  passed  through  a  limiter,  which  removes  any  amplitude  modulation  in 
the  signal  [17].  However,  a  limiter  is  not  needed,  when  k^  is  kept  relatively  small.  For 
the  purposes  of  sub-band  modulation,  the  modulation  index  is  also  kept  relatively  small, 
so  the  modulation  bandwidth  can  be  limited  to  the  desired  frequency  range.  A  further 
discussion  of  spectrum  considerations  for  sub-band  modulation  is  presented  in  section 
4.1. 


3.3  Demodulation 

The  received  echo  is  demodulated  to  recover  the  sub-band  modulating  signals, 
before  the  individual  waveform  components  are  processed.  The  demodulation  process 
is  first  analyzed  without  the  inclusion  of  noise,  and  then  followed  by  an  investigation 
of  the  sub-band  demodulation  process  in  the  presence  of  noise.  In  the  same  vein,  the 
demodulation  process  will  first  be  introduced  without  the  inclusion  of  time  delay  and 
Doppler,  and  then  further  analyzed  with  the  effects  of  time  compression  and  delay  taken 
into  account.  This  incremental  breakdown  and  examination  allows  for  all  of  the  aspects 
of  the  received  signal  affecting  the  demodulation  process  to  be  individually  analyzed. 


3.3.1  Demodulation  without  Interference 


The  point  scatterer  echo  signal  component  of  the  received  data  (2.9)  is 


N  \ 

m=l  / 


(3.26) 


yd  =  1  (  or  (^  =  0  )  and  t  =  0.  A  demodulation  process  is  implemented  to  ex' 


tract  the  modulating  sonar  signals 


from  the  carrier  signal,  and 


an  ideal  discriminator  was  implemented  for  this  thesis.  According  to  Van  Trees  [34] ,  “An 
ideal  discriminator  is  simply  a  device  whose  voltage  output  is  proportional  to  the  instan¬ 
taneous  frequency  of  its  input.”  FM  and  PM  demodulation  using  a  discriminator  (or 
discrimination)  is  the  separation  of  the  modulating  signals  from  the  carrier.  Therefore, 
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the  output  of  a  sub-band  modulation  discriminator  is 

N 

=  V<‘>  =  Vp  E 

m=l 

where  is  a  discriminator  constant  [2,  41]. 

The  ideal  discriminator  used  in  this  study  is  based  upon  the  demodulator  imple¬ 
mented  in  the  Matlab  Signal  Processing  Toolbox  [14].  The  use  of  this  demodulator  is 
discussed  in  section  4.3,  but  the  theory  behind  it  is  presented  here.  First,  the  received 
echo  signal  is  made  analytic  with  a  Hilbert  transform,  and  yields: 


Sj(t)  =  s{t)+jSjj{t) 

N  \  /  N 

771=1  /  \  771=1 

=  ^s^a{t)  cos  (wpi  +  V’(«))  +  J  (0  sin 

=  exp  -I-  jV’(i)) 

"  exp(jWgf)  (3.28) 


The  carrier  is  then  extracted  by  multiplying  Sj(t)  by  the  complex  conjugate  of  the  carrier 
signal, 

S2it)  =  Si(<)exp(-jWQt) 

=  exp  exp 

=  y^S/'(t)  (3.29) 


and  it  becomes  a  baseband  signal. 
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Next,  the  sonax  signals  are  retrieved  from  s^{t)  by  taking  the  imaginary  part  of 
it’s  natural  log: 

s^{t)  =  3{ln[s2(t)]} 

expO'^)]  [o(t)exp(iV’(0)])} 

=  Q|ln(ye^|5^ja(t))  +ln(exp[;'  +  9)])^ 

—  'ip{t)+9 

N 

=  k  y  f  it) +  9.  (3.30) 

p  ^  •'m' 
m=l 

The  derivation  of  (3.30)  shows  that  the  effects  of  the  amplitude  modulation  a{t)  are 
discarded  by  the  discriminator.  The  benefits  of  this  are  discussed  in  section  3.3.2.  (3.30) 
retains  the  uniform  phase  component  of  the  point  scatterer,  but  it  is  turned  into  a 
constant  by  the  discrimination  process,  and  can  be  easily  removed  by  filtering.  The 
highpass  filtered  version  of  (3.30)  is  represented  by  Finally,  the  signal  is  multiplied 

^ ,  and  the  discriminator  yields  the  demodulated 

sonar  signals: 

P 

N 

m=l 

which  is  the  same  form  as  (3.27).  The  waveforms  were  intentionally  designed  to  be 
spectrally  separable  using  bandpass  filtering  at  the  output  of  the  ideal  discriminator.  A 
further  discussion  of  the  waveform  design  criteria  is  presented  in  section  4.2. 

When  #  1  and  t  ^  0,  the  analytic  form  of  the  wideband  Doppler  model  of  the 
echo  signal  component  of  the  data  is 


by  a  discriminator  constant  ^ 


s^it)  =  (/3  [t  -  t])  exp  [t  -  r])  exp  {j2TT<f)t). 


(3.31) 
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First,  the  carrier  is  extracted  from  (3.31)  by  multiplying  by  the  complex  conjugate  of 
the  carrier  signal,  as  was  done  in  (3.29). 

S2{i)  =  [siit)]  exp 

=  {0  [<  -  t])  exp  (juj^  [t  -  r])  exp  (;27r<^t)]  exp 

exp  exp  ij2ir<f>t)^  (3.32) 

★ 

(3.32)  illustrates  that  this  process  does  not  entirely  baseband  the  modulating  sig¬ 
nals  when  P  ^  I,  because  the  Doppler  component  of  the  carrier  (*)  is  not  removed. 
This  mismatch  between  the  transmit  carrier  frequency  used  for  demodulation  and  the 
received  carrier  frequency  effected  by  ★  degrades  the  discriminator’s  performance.  *  can 
be  removed  from  the  received  signal  by  demodulating  with  the  complex  conjugate  of  the 
Doppler  shifted  carrier.  However,  this  assumes  that  the  Doppler  shifted  carrier  frequency 
is  known.  This  is  not  a  practical  assumption  for  a  demodulator  in  a  sonar  signal  pro¬ 
cessing  system,  since  the  Doppler  shift  should  be  estimated  from  the  received  signal  and 
is  not  known  prior  to  signal  reception.  Demodulation  with  a  mismatched  and  a  matched 
carrier  are  investigated  in  the  experiment  and  analysis  section  of  this  thesis  (section  4.4 
cind  5.1).  The  practical  advantages  cind  disadvantages  of  carrier  frequency  matching  in 
demodulation  or  not  will  be  further  discussed  in  these  sections.  However,  it  will  be  shown, 
that  although  a  carrier  mismatch  in  the  demodulation  process  with  an  ideal  discrimina¬ 
tor  does  result  in  decreased  demodulator  performance,  it  does  not  severely  degrade  the 
performance.  Furthermore,  the  resulting  performauce  remains  constant  even  as  the  dif¬ 
ference  between  the  demodulation  frequency  and  the  received  carrier  frequency  increases. 
Therefore,  the  ideal  discriminator  is  still  a  reasonably  effective  sub-band  demodulator, 
even  when  the  carrier  frequency  is  mismatched.  Since  the  ideal  discrimination  process 
is  identical  except  for  the  carrier  frequencies  in  the  cases  of  matched  and  mismatched 
carriers,  the  demodulation  process  will  be  described  here  taking  both  processes  into  con¬ 
sideration.  Following  this  discussion,  a  method  of  coupling  the  demodulation  process 
with  hypothesis  testing  in  order  to  implement  frequency  matching  is  presented.  The 
further  comparison  between  the  different  demodulation  processes  is  left  to  the  analysis 
of  this  thesis  (section  5.1). 
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Equation  3.32  with  the  mismatched  carrier  is  expanded  before  the  demodulation 
process  continues,  so  s^it)  becomes 

=  \/^[|  ‘S’J  exp  (ji9)]  [a  [t  -  r])  exp  (j  [tp  (/?  [f  -  t])])]  exp  exp  O'27r0t) 

=  y^|5ja(/3[t-r])exp(;  [^(/?[t-Tl) +0-a;^r  +  27r</.tJj.  (3.33) 

The  imaginary  paxt  of  the  natural  logarithm  of  s^{t)  is  then  separated  from  the  rest  of 
the  signal,  which  yields  the  modulating  signals,  as  shown  by  (3.30). 


=  ^>{ln  |5^|  a{^[t-  t])]  +  In  [exp  {j  {p  [t-T])+9-  uj  +  27r<^t])]  J 

=  ‘\l}{P\t  -  t])  +  9  —  u^T +  2n<f)t 
N 

=  KYI  [*  -  '^D  +  ^  -  +  27r<^t  (3.34) 

m=l  ' - - - ' 

* 

(3.34)  contains  the  modulating  signals,  along  with  a  linecu:  term  and  two  constants 
(*)  due  to  the  time  dilated  and  delayed  carrier  and  the  phase  of  the  point  scatterer. 

9  and  do  not  have  a  significant  effect  on  the  demodulated  signal,  and  some  of  the 
effects  o{2'K(f>t  can  be  removed  with  a  highpass  filter.  When  the  filter  is  applied  to  (3.34), 
it  yields  essentially  the  modulating  signals  , 

771=1 

The  modulation  index  is  then  removed  from  the  signals  to  yield  the  output  of  the  dis¬ 
criminator, 

N 

(3.35) 

m=l 


the  carrier  frequency  happened  to  be  precisely  matched  by  the  demodulator,  the  highpass 
filter  would  yield  only  the  time  delayed  and  dilated  modulating  signals  without  any  degradation. 
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These  signals  can  then  be  separated  from  each  other  using  band-pass  filters  and  made 
analytic  for  further  processing.  Since  a  broad  bandwidth  is  one  of  the  characteristics 
of  sub-band  modulation,  the  wideband  Doppler  model  is  assumed  for  the  echo  signal 
before  demodulation.  However,  the  narrowband  assumption  may  be  implemented  in  the 
further  processing  of  some  of  the  demodulated  sonar  signals,  if  they  meet  the  criteria  set 
in  (2.14). 

Before  the  sonar  signals  are  processed  according  to  the  methods  discussed  in 
section  2.3,  they  need  to  be  frequency  translated.  Although  the  signals  have  been  de¬ 
modulated  to  their  premodulation  frequencies,  they  have  experienced  time  dilation  at  a 
frequency  near  the  carrier  frequency.  This  results  in  a  Doppler  shift  that  is  much  larger 
than  if  the  signal  had  not  been  sub-band  modulated  at  all.  This  frequency  shift  becomes 
evident  by  comparing  only  one  of  the  sub-band  modulated/demodulated  sonar  signals 
l)  same  signal  transmitted  without  the  use  of  sub-band 

modulation  ^  center  frequency  of  and  experiences  a  Doppler 

shift  of 

2u,r7, 


Although  has  the  same  center  frequency,  it  was  phase  modulating  a  carrier  at  77^, 

where  <C  rj^.  Therefore,  the  Doppler  frequency  of  ^(t)  after  demodulation  is 


2v,tj 
I  'o 

c 


(3.37) 


even  though  the  center  frequency  of  the  signal  is  rjy  This  presents  a  problem,  when 
processing  ^(t)  with  a  matched  filter,  because  the  Doppler  frequency  of  the  demodu¬ 
lated  signal  is  no  longer  representative  of  the  Doppler  velocity.  The  estimated  scatterer 
velocity  is  instead  of  v^.  Therefore,  the  demodulated  and  separated  signals  are 
shifted  in  frequency  up  to  the  carrier  frequency  to  account  for  this  processing  anomaly. 
This  is  achieved  by  multiplying  the  signal  by  the  carrier  signal.  This  shift  in  frequency 
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is  known  as  a  heterodyne  shift,  and  the  demodulated  signal  becomes 


exp 

=  ipt)  exp  (j  +  Wq]  t)  exp  ^  tj ,  (3.38) 

where  +  Wq  w  Wq  because  o/^  ^  ^0’  Doppler  shift  in  the  demodulated  signal 

is  representative  of  the  velocity  of  the  scatterer.  When  processed  with  a  matched  filter 
bank  of  time  dilated  and  delayed  replicas  centered  ciround  the  correct  Doppler 

Q 

values  are  estimated  .  It  will  be  shown  in  section  4.2,  that  this  frequency  shift  in  the 
demodulation  process  is  taken  into  consideration  in  the  waveform  design  process. 

Ironically,  the  ideal  discriminator  is  not  always  the  best  choice  for  a  demodulator. 
Prior  knowledge  of  the  carrier  frequency  and  deviation  constant  zire  required  for  effective 
demodulation.  Other  demodulators  can  be  more  advantageous  in  certain  situations,  such 
as  the  phase-lock  loop  [30].  The  phase- lock  loop  is  a  feedback  demodulator,  that  im¬ 
plements  a  voltage-controlled  oscillator  to  track  the  frequency  deviations  of  the  received 
signal.  A  voltage  proportional  to  the  frequency  deviation  between  the  received  signal 
and  a  receiver  generated  sinewave  are  applied  to  the  voltage-controlled  oscillator.  This 
shifts  its  frequency  deviations  to  approximate  the  received  signal  [27,  12,  28,  34,  41].  The 
frequency  deviation  of  the  received  signal  is  measured  by  implementing  this  feedback  de¬ 
modulator,  and  the  phase-lock  loop  demodulates  a  received  FM  or  PM  signal  without 
prior  knowledge  of  the  carrier  frequency.  The  phase-lock  loop  also  exhibits  good  noise 
performance,  because  the  frequency  deviations  are  constantly  tracked  [28].  However, 
phase-lock  loop  (and  other  demodulators)  require  a  much  more  complex  and  rigorous 
implementation.  Since  the  object  of  this  research  is  to  investigate  the  application  of 
phase  modulation  in  sonar,  the  investigation  of  more  complex  demodulation  schemes 
will  be  left  for  future  research. 

One  method  of  coupling  the  ideal  discriminator  with  frequency  matching  as  part 
of  the  hypothesis  testing  process  is  explored  in  this  thesis.  Instead  of  demodulating  the 
received  signal  with  only  one  carrier  frequency  and  then  matched  filtering,  the  received 

Only  is  used  in  (3.38)  instead  of  /^(/?[t  —  t]),  because  the  effects  of  time  delay  were 

removed  to  simplify  processing.  This  was  done  using  methods,  that  will  be  presented  in  section 
4.3. 


dh 


i(*)  = 


f^Wt)  exp 


(■'“lO 


exp  j 
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signal  is  demodulated  with  a  time  dilated  carrier,  whose  time  dilation  correspond  with 
the  time  dilation  of  the  replica,  against  which  it  would  be  correlated  in  the  matched 
filter  process.  Therefore,  (3.32)  becomes 

S2W  =  [siW]exp(-ja;gt) 

=►  [sj(t)]  exp 

=  [s^(t)]  exp  exp  .  (3.39) 

This  results  in  an  array  of  demodulated  signals,  each  demodulated  with  a  different  carrier 
exp  This  means,  that  the  demodulated  sigiuils  before  spectral  separation 

(from  (3.34))  become 

N 

s^{t)  =  ^  Y1  fm  +  27r(^t  -  2'K^t.  (3.40) 

m=l 

Therefore,  the  carrier  time  dilation  is  canceled  out  when  ^  which  occurs  when  the 
received  carrier  frequency  is  matched  by  the  discriminator.  This  array  of  demodulated 
signals  could  then  each  be  spectrally  separated,  and  the  resulting  arrays  of  sub-band 
signals  could  then  be  matched  filtered. 

3.3.2  Demodulation  in  the  Presence  of  Interference 

The  demodulation  process  becomes  even  more  complex,  when  interference  is  also 
tciken  into  consideration.  The  echo  signal  embedded  in  additive  interference  is  repre¬ 
sented  by 


r(t)  =  s{t)  -I-  n(t) 

=  ^Kl  cos  +  0  +  V’(t)j  +  ■ 


(3.41) 


The  interference  term  (  n{t) )  can  be  either  reverberation  or  noise,  and  the  demodulation 
process  is  identical  in  the  presence  of  either  one.  n{t)  in  this  discussion  is  additive  white 
Gaussian  noise,  where  a  (t)  is  the  Rayleigh  distributed  noise  envelope  and  il)  {t)  is  the 
uniformly  distributed  phase  of  the  noise  [41].  An  effective  way  to  illustrate  the  impact 
of  noise  on  a  sub-band  modulating  signal  is  to  represent  the  received  signal  in  vector 
notation.  The  vector  of  magnitude  rotates  at  the  carrier  frequency  in  a 

noiseless  environment,  and  that  rotation  is  phase  modulated  by  the  message  signals. 
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d  is  the  phase  shift  imposed  upon  the  echo  signal  by  the  scatterer  (2.7),  and  it  is 
set  to  zero  here  to  simplify  the  vector  addition.  The  addition  of  interference  results  in 
the  vector  addition  between  n{t)  and  the  echo  signal  at  the  receiver.  This  yields  a  new 
vector  amplitude  and  phase,  as  seen  in  figure  3.7. 

The  received  signal(3.41)  is  rewritten  with  vector  addition  as 


II 

s 

s 

COS 

Wpt -1-^(0 j  -l-a^(t)cos 

-1-  i>{t)  + 

= 

S 

s 

cos 

+  V'(0]  +a^(t)cos 

+  rp{t)^  cos  -  ■^('(t)] 

—  a 

^{t)  sin  -1-  V’(0|  sin 

V'„(<)  -  i’it)] 

=  1-5^1  +  a„(<)  cos  [V’^(t)  -  cos 

-a^{t)  sin  -'?/'(*)]  sin  [wpt  +  V’(<)] 

=  R{t)  cos  (3.42) 

,  where 


(3.43) 


and 


'ip{t)  +  tan 


(  o„(<)sin 

] 

Vv/% 

5 

s 

-h  a  (t)  cos 
n'  ' 

- 

-m]) 

(3.44) 


[2,  9,  27,  28,  41].  This  notation  illustrates  the  signal  resulting  from  the  addition  of 
interference  and  the  received  echo,  and  this  combined  signal  is  the  one,  that  gets  demod¬ 
ulated. 

The  effectiveness  of  the  sub-band  demodulation  process  depends  solely  upon 

(3.44),  since  the  discriminator  disregards  (3.43),  as  seen  in  section  3.3.1.  If  a  high 
signal-to-interference  ratio  at  the  receiver  (SIR^)  is  present  (ojj(^)  <  \/%  1*^31)’ 
the  angular  variations  in  the  noise  will  result  in  minimal  effects  on  the  angle  of  the 
message  signal,  as  seen  in  figure  3.7.  Therefore,  (3.44)  Ccin  be  rewritten  as 


a  (t) 


’’Pjlii)  =  t{t)  + 


sin  [i’Ji)  -  i>{t)  , 


(3.45) 
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Fig.  3.7.  Vector  Addition  of  Echo  Signad  ^s(t)  =  cos  +  V’(*)])  and  Noise 

^n(t)  =  cos  +  '^^(^)])  at  a  Moment  in  Time  (High  SIR^) 


where  ^  — r  is  small  compared  to  and  the  phase  of  the  received  signal  is  domi- 
nated  by  the  message  signed  [41]. 

However,  the  received  signad  dramatically  changes  as  the  signal- to- interference 
ratio  becomes  smadl.  The  noise  signal  becomes  the  dominamt  part  of  the  received  signal 
and  the  angular  vairiations  of  the  echo  signal  will  have  only  a  very  small  effect  on  its 
angle  for  a  low  signed- to-interference  ratio  at  the  receiver  )>  ^ 

in  figure  3.8  [9].  Therefore,  (3.42)  is  rederived  with  the  noise  signal  as  the  dominating 
signad,  and  (3.44)  is  rewritten  as 

^^(t)  =  sin  [^^(t)  -  V’(<)] ,  (3.46) 


^/S-S 
Vo  s 


where  ^  is  small  compared  to  Tp  {t)  [41]. 

Although  this  seems  like  a  drastic  difference  in  demodulator  performance,  an¬ 
gle  modulation  (and  sub-band  modulation)  offers  significant  demodulation  improvement 
over  amplitude  modulation  in  a  noisy  environment.  The  magnitude  of  the  received  signal 


48 


Fig.  3.8.  Vector  Addition  of  Echo  Signal  ^s(t)  =  cos  and  Noise 

^n(t)  =  a^{t)co9  at  a  Moment  in  Time  (Low  SIR^) 


varies  linearly  with  the  signal-to-interference  ratio  [9].  Figure  3.9  shows  the  lineair  in¬ 
crease  in  magnitude  of  the  received  signal  as  compared  to  the  magnitude  of  only  the  echo 
signal.  The  amplitude  difference  increases  in  inverse  proportion  to  SIR^,  which  in  turn 
adversely  2iffects  demodulation  of  an  amplitude  modulating  echo  signal.  Although,  this 
strongly  effects  the  performance  of  AM,  it  is  not  detrimental  to  sub-band  modulation, 
since  the  PM  demodulation  disregards  the  amplitude  of  the  received  signal.  This  is  com¬ 
parable  to  the  reception  of  an  AM  radio  signal,  which  does  not  fade  out  right  away  when 
the  SIR^  gets  small,  but  instead  gets  worse  and  worse  as  the  SIR^  decreases.  However, 
interference  affects  the  phase  of  the  received  signal  differently.  Interference  has  very  lit¬ 
tle  effect  on  the  phase  of  the  received  signal  until  the  SIR^  gets  small.  When  the  SIR^ 
gets  sufficiently  small,  the  noise  quickly  dominates  the  phase  component  of  the  received 
signed,  which  is  known  as  the  threshold  effect  [29].  Figme  3.10  shows  the  drastic  change 
in  phase  with  decreasing  SIR.  as  compared  to  the  phase  of  only  the  echo  signal.  This 
primarily  affects  angle  modulation,  because  it  depends  so  heavily  upon  the  phase  of  the 
received  signal  for  effective  demodulation.  This  drastic  change  in  phase  with  decreasing 
SIR^  is  comparable  to  FM  radio.  FM  radio  has  a  much  clearer  reception  than  AM  over 
a  broeider  SIR.  range,  but  then  quickly  drops  out,  when  the  SIR.  has  decreased  below 
a  certain  point. 
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Fig.  3.9.  Maximum  Amplitude  DifiFerence  between  the  Received 
Signal  (Signal  plus  Noise,  R{t))  and  the  Echo  Signal  (|s(t)|) 
with  Respect  to  Received  Signal-to-Interference  Ratio  (Constant 
Noise  Level  as  Signal  Level  Increases) 


Fig.  3.10.  Maximum  Phase  Difference  between  the  Received 
Signal  (Signal  plus  Noise,  and  the  Echo  Signal  (^(t)) 

with  Respect  to  Received  Signal-to-Interference  Ratio  (Constant 
Noise  Level  as  Signal  Level  Increases) 
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Chapter  4 


Experiment 


The  experiment  portion  of  this  thesis  is  divided  into  five  section  to  better  present 
the  implementation  of  sub-band  modulation  techniques  in  active  sonar.  First,  the  nu¬ 
merous  factors  contributing  to  waveform  design  are  discussed.  Second,  the  specifications 
of  the  waveforms  used  in  this  thesis  are  presented,  along  with  the  motivation  behind 
the  different  combinations  of  sonax  waveforms.  Third,  the  implementation  of  sub-band 
modulation  process  in  Matlab  is  summarized.  Fourth,  ein  overview  the  simulations  used 
in  this  thesis  is  given,  paying  particular  attention  to  the  setup,  structure  and  goals  of 
the  simulation  experiments.  Finally,  the  conditions  and  goals  behind  the  in-water  ex¬ 
periment  section  of  this  thesis  are  discussed. 

4.1  Waveform  Design  Considerations 

Waveform  design  is  an  important  part  of  sub-band  modulation.  Many  spectrum 
considerations  need  to  be  taken  into  account  when  designing  a  sub-band  modulating 
waveform,  because  of  the  numerous  limitations  imposed  by  the  transmitting  transducer, 
the  demodulation  process  and  the  effects  of  time  compression.  Therefore,  the  factors 
that  go  into  waveform  design  are  expanded  upon  here. 

The  transmission  of  a  sub-band  modulated  waveform  is  achieved  with  a  trans¬ 
mitting  transducer  (or  array  of  transducers),  and  the  signal  is  affected  by  the  frequency 
response  of  this  transmitter.  Each  sonar  transmitter  has  a  set  operating  bandwidth,  and 
the  spectnun  of  the  sub-band  modulating  waveform  needs  to  fall  within  this  bandwidth 
to  ensure  proper  transmission  into  the  underwater  environment.  Carson’s  rule  (3.23) 
should  be  implemented  when  designing  the  sonar  waveform  and  selecting  the  deviation 
constant,  so  the  spectrum  of  the  sub-band  modulating  waveform  will  fall  within  the 
transmitter  ba,ndwidth. 

The  effects  of  time  dilation  also  need  to  be  taken  into  consideration.  The  trans¬ 
mitted  signal  is  frequency  shifted  by  the  Doppler  effect  (2.13),  and  if  this  shift  is  large 
enough,  the  message  signals  can  be  shifted  out  of  the  frequency  range  of  the  receive 
transducers.  When  designing  a  sub-band  modulating  waveform,  precautions  need  to  be 


51 


taken  in  order  to  ensure  that  the  signal  will  not  be  Doppler  shifted  out  of  the  receive 

bandwidth  B..  The  deviation  constant  k  can  be  selected  to  meet  the  desired  band- 
t  p 

width  criteria  by  making  general  assumptions  about  the  maximum,  closing  line-of-sight 
velocity  of  potential  scatterers  (  implementing  Carson’s  rule  (3.23), 

B. 

k  <  7: - - - 1,  (4.1) 

P  2a;  '  ^ 

^Imax 

where  is  the  bandwidth  of  the  transducer  and 


j  —  (jj^  I  1  + 

^Imax  0  V 


Therefore,  the  bandwidth  of  the  modulating  waveform  can  be  designed  so  it  still  falls 
within  the  receive  bandwidth  even  in  the  case  of  maximum  time  dilation  or  compres¬ 
sion.  Typically,  the  receive  transducer  has  a  bandwidth  wide  enough  to  account  for  the 
maximum  frequency  shift,  experienced  by  signals  within  the  transmitter’s  bandwidth. 
However,  the  extra  design  criteria  are  taken  in  account  as  an  additional  precaution. 

The  transmitter  response  is  not  always  flat  over  it’s  entire  bandwidth.  Either 
only  the  flattest  section  of  the  frequency  spectrum  of  the  transducer  is  utilized,  or  the 
deviation  constant  and  the  frequencies  of  the  sonar  waveforms  must  be  selected  in  such 
a  way  that  the  sub-band  information  does  not  fall  within  the  dips  of  the  transmitter’s 
response.  If  portions  of  the  modulating  sonar  signals  fall  in  dips  in  the  transmitter’s 
frequency  response,  then  the  information  at  those  frequencies  is  broadcast  at  a  lower 
level  than  the  rest  of  the  signal,  which  may  present  problems  for  demodulation. 

The  selection  of  the  specific  sonar  waveforms  to  be  used  in  sub-band  modulation 
is  an  involved  process.  Besides  choosing  a  combination  of  waveforms  that  can  estimate 
the  desired  sonar  parameters  (as  discussed  in  sections  2.3  and  2.4),  the  waveforms  are 
designed,  so  they  can  be  spectrally  separated  after  demodulation.  When  the  sonar 
waveforms  are  combined  linearly  as  the  phase  component  of  a  csirrier  (3.16),  they  must 
be  spaced  far  enough  apart  in  frequency,  that  they  can  be  separated  using  filtering 
techniques.  The  filters  used  to  separate  these  signals  also  need  to  have  wide  enough 
passbands  to  allow  for  the  potential  Doppler  shift  of  the  sonar  signals  (2.13), 

where 

^filter  m  -  m  1  c 


(4.3) 
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Therefore,  waveforms  adjacent  in  frequency  should  be  spectrally  spaced  according  to 
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(3.20)  shows,  that  the  the  sub-bands  of  the  modulation  lie  at  the  frequencies 
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The  center  frequencies  of  the  sonar  waveforms  need  to  be  chosen,  so  there  are  no 

significant  sidebands  at  frequencies  within  the  bandwidths  of  the  filters  used  for  post¬ 
demodulation  signal  separation.  Prominent  sidebands  may  appear  as  additional  peaks 
in  the  matched  filter  response  and  should  be  avoided  to  ensure  effective  demodulation. 

Filtering  is  used  in  the  demodulation  process  to  remove  unwanted  signal  compo¬ 
nents  from  the  sonar  signals  (sections  3.3.1  and  3.3.2),  as  well  as  to  spectrally  separate 
the  multiple  sonar  signals.  The  sonar  signal  with  the  lowest  center  frequency 
should  be  designed,  so  that  it’s  frequency  is  far  enough  away  from  these  components 
to  allow  for  their  removal.  This  filter  is  a  very  low  frequency  bandstop  filter  with  a 
bandwidth  of  The  ctilmination  of  the  constraints  on  the  sonar  waveforms’  spectra 
implemented  in  sub-band  modulation  is  illustrated  in  figure  4.1. 


4.2  Waveform  Specifications 

Severed  different  combinations  of  modulating  waveforms  were  implemented  through 
out  the  course  of  this  research  to  verify  the  effectiveness  of  sub-band  modulation.  Specifi¬ 
cally,  one  waveform  was  used  through  out  the  majority  of  the  research,  and  several  other 
waveforms  were  used  in  parts  of  the  research  to  verify  that  the  sub-band  modulation 
performance  was  independent  of  the  types  of  sonar  waveforms  used, 
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frequency  (radians/sec.) 


Fig.  4.1.  Spectral  Constraints  of  Waveform  Design:  /^(<)  Waveform  “m”  with  oj 


m 


Center  Frequency,  ^  Filter  for  Separating  Waveform  “m”  from  Other  Modulat¬ 

ing  Waveforms,  B  Bandwidth  of  Modulation,  Bandwidth  of  Receiver  and  B^  Low 
Frequency  Bandstop  Filter 


The  primary  waveform  used  in  this  research  is  a  combination  of  a  CW  and  a 

BPSK, 

where  the  center  frequency  of  the  BPSK  is  with  a  mainlobe  bandwidth  of 

(see  figure  4.2a).  w,  and  u;_  are  the  basebanded  center  frequencies.  The  main  motiva- 
tion  behind  this  waveform  combination  is  to  simultaneously  implement  a  waveform  with 
high  frequency  and  time  resolution  along  with  waveform  with  poor  time  resolution.  This 
combination  uses  the  CW  to  make  a  detection  estimation  and  to  make  generalizations 
about  scatterer  parameters.  Then,  the  desired  high  resolution  scatterer  parameters  are 
extracted  from  the  matched  filter  response  of  the  BPSK  using  a  more  focused  bank  of 
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replicas.  This  waveform  combination  allows  for  a  coarse  approximation  of  the  echo  pa¬ 
rameters  made  with  a  low  resolution  signal  to  ease  the  computational  intensity  associated 
with  the  high  resolution  signal,  as  discussed  in  section  2.4. 

This  waveform  combination  was  also  chosen  to  verify  and  investigate  the  band¬ 
width  constraints  of  sub-band  modulation.  Since  the  BPSK  has  a  reasonably  large  band¬ 
width,  it’s  sub-band  modulation  with  any  other  sonar  waveform  produces  a  wideband 
waveform.  Therefore,  the  BPSK  waveform  is  bandpass  filtered  to  keep  the  modulating 
bandwidth  within  the  transducer’s  bandwidth  constraints.  It  is  found,  that  95.45%  of  the 
signal  energy  in  this  BPSK  is  contained  within  the  mainlobe  and  the  first  two  sidelobes 
of  the  sine  function.  The  overall  bandwidth  of  the  BPSK  can  be  reasonably  reduced 
without  compromising  the  waveform  energy  by  bandpass  filtering  out  some  sidelobes. 
This  also  allows  for  wider  spectral  separation  between  the  two  waveforms,  in  order  to 
facilitate  easier  signal  separation  in  the  demodulation  process. 

Several  other  signal  combinations  are  used  in  parts  of  this  research,  to  verify 
that  the  sub-band  modulation  performance  is  independent  of  sonar  signal  type.  One 
waveform  combination  is  the  CW  and  LFM 


(4.6) 


where  is  a  LFM  centered  at  with  a  bandwidth  of  as  seen  in  figure  4.2b. 

The  combination  of  a  signal  with  good  low  Doppler  performance  with  a  signal  with  good 
high  Doppler  performance  is  implemented  to  verify,  if  sub-band  modulating  signals  could 
yield  good  Doppler  performance  independent  of  scatterer  velocity.  Another  waveform 
combination  is  four  CWs  of  different  frequencies, 


(4.7) 


as  seen  in  figure  4.2c.  This  waveform  combination  is  used  to  investigate  the  effects  of 
the  number  of  sub-band  waveforms  on  the  modulation  process.  Finally,  a  combination 
of  2  sinewaves  of  duration  T  and  two  sinewaves  of  duration  ^  is  investigated.  This 


[d]  norm.  dB  [c]  norm.  dB  [b]  norm.  dB  [a]  norm.  dB 
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is  seen  in  figure  4.2d.  This  waveform  combination  investigates  the  effects  of  using  wave¬ 
forms  with  different  durations  on  sub-band  modulation. 

Realistically,  all  of  the  waveforms  used  in  this  thesis  would  have  center  frequencies 
of  a;  -I-  a;  instead  of  just  cj  ,  if  they  were  transmit  without  the  use  of  sub-band 
modulation.  The  waveforms  are  basebanded,  before  they  are  sub-band  modulated,  so 
the  waveforms  in  the  first  sub-band  of  the  carrier  would  be  at  the  same  center  frequency 
as  those  transmitted  without  sub-band  modulation.  This  was  done  to  minimize  the 
bandwidth  of  the  modulation  and  to  ease  processing  complexity,  by  having  the  sub¬ 
band  modulated  and  the  equivalent  unmodulated  signals  experience  essentially  the  same 
Doppler  shift. 

Once  the  sonar  waveforms  are  selected,  the  carrier  frequency  and  deviation  con¬ 
stant  of  the  sub-band  modulation  cire  carefully  considered.  The  most  stringent  design 
parameter  in  this  research  is  the  bandwidth  of  the  transmitting  transducer  implemented 
in  section  4.5.  Since  is  to  be  the  only  waveform  combination  implemented  in  all 
parts  of  the  experimental  stage  of  this  research,  it’s  design  is  highly  influenced  by  these 
bandwidth  constraints.  A  large,  relatively  flat  portion  of  the  transducer’s  frequency  re¬ 
sponse  is  utilized  in  this  research,  so  the  carrier  frequency  and  modulation  bandwidth 
are  selected  to  fit  within  this  frequency  range.  The  carrier  frequency  Wq  is  selected  to  fall 
in  the  center  of  the  desired  transmit  band,  and  Carson’s  rule  determined  the  deviation 
constant,  that  will  keep  the  maximum  frequency  shift  due  to  the  scatterer  velocity  within 
the  transmission  bandwidth  (4.1).  The  value  of  the  modulation  index  =  0.5^  is  also 
selected  so  that  the  first  set  of  sub-bands  are  much  more  prominent  them  most  of  the 
other  sub-bands.  Therefore,  the  sub-band  modulating  signals  are  essentially  represented 
by  the  first  set  of  sub-bands.  Although,  these  first  sub-bands  have  a  low  magnitude  as 
compared  to  the  carrier,  is  kept  small  ^k^  <  f )»  because  a  higher  modulation  index 
drastically  increases  the  modulation  bandwidth  along  with  the  sub-band  magnitudes, 
as  discussed  in  section  3.2.1.  Since  the  bandwidth  of  the  modulation  is  already  wide, 
due  to  the  inclusion  of  a  BPSK,  the  modulation  index  is  kept  relatively  small  to  avoid 
sub-band  modulation  outside  of  the  transmission  bandwidth.  The  frequency  spectrum 
of  this  sub-band  modulated  signal  can  be  seen  in  figure  4.3. 
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Filtered  and  Sub-band  Modulated  Waveform  \j/  (t)  before  Transmission 
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Fig.  4.3.  Frequency  Spectrum  of  the  Sub-band  Modulated  BPSK  and  CW  Waveforms 
in  V’g(t)  with  Modulation  Index  —  0.5,  =  0.18^  and  B  =  0.06 Immediately 

before  Transmission  (Normalized  Frequency  to 
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Although  was  chosen  according  to  the  design  criteria  in  section  4.2  for  use  with 
signal  combination  it  was  subsequently  used  for  all  four  test  signal  combinations, 

because  they  had  very  similar  bandwidths.  This  modulation  index’s  effectiveness  with 
all  four  test  signal  combinations  was  confirmed  through  a  series  of  simulations,  where 
different  deviation  constants  were  compared.  One  possibility  for  future  work  may  be  to 
investigate  a  sub-band  modulating  signal  with  a  narrower  modulation  bandwidth,  so  a 
more  extensive  analysis  of  the  modulation  index  may  be  made. 

4.3  Matlab  Implementation 

This  section  of  the  thesis  is  devoted  to  the  application  of  sub-band  modulation 
theory  as  a  series  of  signal  processing  programs.  The  waveform  generation,  sub-band 
modulation,  demodulation,  signal  manipulation  and  echo  processing  in  this  research  are 
all  implemented  in  Matlab  version  5.3  with  Signal  Processing  Toolbox  [14].  This  section 
provides  a  general  overview  of  the  sub-band  modulation  processes  in  Matlab,  which  will 
then  be  applied  in  the  simulation  and  in-water  experiment  sections  of  this  thesis  (sections 

4.4  and  4.5). 

Most  of  the  waveforms  in  this  research  are  generated  using  standard  mathematical 
functions  in  Matlab.  The  one  exception  is  the  BPSK,  which  is  generated  with  a  function 
created  by  Anthony  J.  Cutezo  [6].  These  discrete  waveforms  are  generated  at  a  sampling 
rate  of  tj  ,  where  uj  w-  w  • 

s  3  [J  m 

Amplitude,  phase/sub-band  and  frequency  modulation  are  achieved  using  the 
“modulateO”  command  from  the  Signal  Processing  Toolbox.  The  source  code  for  this 
function  was  compared  to  the  modulation  theory  discussed  in  section  3.2  directly  applied 
in  Matlab,  and  was  found  to  yield  identical  results.  Since  the  “modulate()”  function  is 
less  computationally  intensive  than  the  independently  developed  modulation  functions, 
it  is  used  through  out  this  research.  This  applies  to  the  “demod()”  function  for  de¬ 
modulation,  as  well.  Matlab  also  has  numerous  filtering  functions,  that  can  be  used  to 
separate  the  demodulated  echo  signals  spectrally  or  to  represent  the  bandwidth  of  the 
transducer. 

The  frequency  matching  of  the  carrier  for  demodulation  is  implemented  in  Matlab 
with  a  fast  Fourier  transform  (fft)  of  the  scatterer  return  portion  of  the  received  time 
series.  The  narrowband  frequency  shift  of  the  carrier  signal  can  be  approximated  from 
the  resulting  frequency  spectrum.  The  rough  estimation  of  the  carrier’s  narrowband 
frequency  shift  is  a  close  enough  approximation  of  the  time  dilation  of  the  carrier,  that  it 
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can  be  used  with  the  ideal  discriminator.  However,  this  method  of  frequency  matching  is 
only  effective  in  high  receive  signal-to-interference  situations,  where  the  scatterer  return 
is  discernible  from  the  rest  of  the  received  time  series. 

The  processing  of  the  individual  sonar  signals  is  also  implemented  in  Matlab. 
After  the  data  is  read  into  Matlab,  the  time  series  data  is  displayed  to  verify  that  the 
echo  signal  is  present  in  the  data.  Since  the  overall  data  set  can  be  several  seconds 
long  at  a  high  sampling  frequency,  the  echo  return  is  cropped  from  the  received  data 
for  processing.  Processing  time  and  complexity  are  greatly  alleviated  by  only  processing 
the  echo  return  portion  of  the  received  data.  Then,  the  echo  location  in  the  received 
data  is  noted  and  applied  to  the  matched  filter  responses  to  achieve  the  desired  range- 
to-scatterer  estimations.  This  method  does  not  work  as  well  in  the  case  of  low  signal- 
to-interference  ratio.  When  the  signeil-to-interference  ratio  is  so  low,  that  the  scatterer 
return  could  not  be  discerned  from  the  reverberation,  alternative  methods  were  used  to 
locate  it.  A  priori  knowledge  of  the  scatterer’s  location  was  utilized  in  the  simulations, 
and  observations  recorded  during  data  collection  aided  in  locating  scatterer  returns  in 
the  in-water  data.  This  did  help  in  locating  the  scatterer  returns  in  the  time  series  with 
low  signal-to-interference  ratio  levels.  However,  it  is  not  a  practical  method  to  use  in 
a  “real”  sonar  detection  scenario,  because  it  relies  on  already  knowing  the  scatterer’s 
location,  before  it  is  estimated.  Although,  such  information  is  not  available  in  “real” 
detection  scenarios,  the  purpose  of  this  resecirch  is  to  verify  that  echo  parameters  can  be 
estimated  from  sub-band  modulating  signals,  which  is  still  achieved  using  this  process. 
Therefore,  it  is  reasonable  to  apply  in  this  research.  Topics  for  future  research  may 
include  more  extensive  processing  techniques,  that  allow  for  processing  the  entire  set  of 
received  data  more  effectively  and  efficiently. 

4.4  Simulations 

Two  sets  of  simulations  were  completed  as  part  of  this  research.  One  is  a  simple 
ocean  setup  to  verify  the  entire  sub-band  modulation  process,  and  the  other  is  a  sim¬ 
ulation  of  the  in-water  experiment  to  better  understand  the  field  results.  The  former 
is  discussed  here,  while  the  later  will  be  presented  in  section  4.5  along  with  the  actual 
in-water  portion  of  this  research  for  better  compjirison  between  the  two. 

The  sonar  simulations  in  this  research  were  created  using  methods  outlined  in 
[8],  which  will  be  referred  to  as  the  synthetic  data  generation  process  (or  SDG)  in  the 
remainder  of  this  thesis.  The  SDG  process  simulates  sound  propagation,  the  ocean 


60 


environment,  ocean  surfaces,  scatterers,  scatterer  echos,  reverberation  and  background 
noise.  The  SDG  process  allows  for  easy  implementation  with  Matlab,  so  the  signals  can 
be  generated  and  modulated  in  Matlab,  sent  to  the  SDG  process  for  signal  propagation 
simulation,  and  then  read  back  into  Matlab  for  demodulation  and  processing. 

The  purpose  of  this  portion  of  the  research  is  to  demonstrate  that  sub-band 
modulation  is  a  viable  option  for  active  sonar.  Therefore,  a  simple  direct  path  ocean 
environment  without  boundary  reflections  is  implemented.  A  single  point  scatterer  at 
the  same  depth  as  the  sonar  platform  is  used  in  a  uniformly  deep  ocean  environment, 
which  allows  for  a  single  speed  of  sound  to  be  used  instead  of  an  entire  sound  velocity 
profile  for  the  different  strata  in  the  ocean.  The  transmit  and  receive  beams  on  the  sonar 
transducer  are  omnidirectional,  so  the  effects  of  beam  directionality  does  not  effect  the 
data.  The  simulation  is  designed  to  be  as  simple  as  possible,  so  the  effects  of  rever¬ 
beration,  target  velocity,  and  target  range  on  sub-band  modulation  can  be  investigated 
without  the  scattering  effects  of  boundary  interactions. 

The  mmn  experiment  consists  of  two  hundred  and  forty  simulations  using  the  com¬ 
bination  BPSK  and  CW  signal  Different  scatterer  velocities  are  implemented 

to  test  the  effect  of  time  dilation  on  demodulation  and  detection  with  a  sub-band  mod¬ 
ulated  signal.  The  varying  scatterer  velocity  also  permitted  for  the  evaluation  of  the 
detection  performance  of  the  sub-band  signals  with  and  without  frequency  matching  in 
the  discriminator.  A  smadl  number  of  these  simulations  were  also  used  to  test  the  per¬ 
formance  of  the  coupled  demodulation,  hypothesis  testing  process.  This  experiment  is 
further  subdivided  according  to  the  signal-to-interference  ratio  at  the  receiver  (^SIR^, 
to  investigate  the  effects  of  interference  on  the  demodulation  and  detection  of  a  sub-band 
modulated  signal.  These  SIR.  values  were  achieved  by  varying  the  location  and  the  scat¬ 
tering  strength  of  the  scatterer,  and  all  of  the  simulations  eveduated  fall  within  a  ±1.5dB 
range  of  the  desired  value.  This  reinge  of  SIR^  values  is  needed,  because  both  the  target 
cind  reverberation  were  randomized  by  the  synthetic  data  generation  process  to  provide  a 
more  realistic  simulation.  The  signal-to-interference  ratio  at  the  receiver  ^SIR^^  is  used 

instead  of  the  signal-to-interference  ratio  at  the  matched  filter  output  ^SIR^ ,  because 
it  permits  the  direct  comparison  between  the  performance  of  a  sub-band  modulating 
signal  and  the  performance  of  the  same  unmodulated  signals  in  equivalent  interference. 
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Two  hundred  and  forty  simulations  using  unmodulated  BPSK  and  CW  signals 
were  also  generated  for  comparison  with  the  sub-band  modulated  signals.  These  simula¬ 
tions  had  the  same  parameters  as  the  previous  simulations,  only  the  signals  were  not  sub¬ 
band  modulated.  The  sonar  waveforms  used  ^  the  BPSK  and  ^  the  CW 

were  the  same  as  the  sub-band  modulating  signals  in  the  previous  simulations  /.  (t)  and  fAt), 
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except  they  are  centered  at  and  respectively  instead  of  being  basebanded. 

This  is  done,  so  the  unmodulated  signals  measured  at  the  receiver  can  be  directly  com¬ 
pared  to  the  received  sub-band  modulated  signals,  which  cire  frequency  translated  up 
to  the  carrier  frequency  once  they  axe  demodulated  (as  discussed  in  section  3.3.1).  The 
results  of  this  comparison  axe  presented  and  analyzed  in  chapter  5.1. 

The  other  experimental  sets  with  the  waveform  combinations  i’JA)  S'nd 
were  carried  out  in  the  same  manner  as  except  fewer  simulations  were  completed 

for  each  set.  Each  set  only  consisted  of  seventy-two  simulations  for  the  sub-band  mod¬ 
ulated  signals  cind  seventy-two  additionaJ  simulations  for  the  equivalent,  unmodulated 
signals.  These  experimental  sets  did  not  need  to  be  as  extensive  as  the  first  one,  since 
their  main  purpose  was  to  investigate  the  influence  of  signal  type  and  number  of  modu¬ 
lating  signals  on  detection  with  sub-band  modulated  signals  and  verify  the  performance 
of  the  first  experiment2d  set. 

4.5  In- Water  Data 

A  small  in-water  experiment  was  also  completed  as  part  of  this  research.  The  sub¬ 
band  modulating  signal  combination  of  i’JA)  used  on  twenty-four  separate  instances 
to  interrogate  a  scatterer  with  Doppler.  The  purpose  of  this  portion  of  the  research  was  to 
experimentally  verify  the  detection  performance  observed  in  the  simulations.  Although 
not  enough  interrogations  were  completed  to  generate  detection  statistics,  the  data  was 
used  to  validate  the  results  gathered  in  section  4.4. 

The  ocean  environment  was  similar  to  the  simulation  set  up  discussed  in  section 
4.4,  however  the  sound  propagation  was  much  more  complex.  Since  the  sound  propaga¬ 
tion  could  not  be  limited  to  the  direct  path  between  the  sonar  platform  and  the  scatterer, 
surface  and  bottom  bounces  were  also  included  in  the  received  data.  A  varying  sound 
velocity  profile  also  affected  the  propagating  sonar  signeils,  as  did  the  high  volume  at¬ 
tenuation  in  the  medium.  All  of  these  factors  contributed  in  creating  a  very  strong 
reverberation  level.  Although  the  ocean  environment  was  not  too  different  from  the  one 
used  in  the  simulations,  the  added  complexity  of  the  reed  ocean  environment  and  higher 
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reverberation  levels  decreased  the  detection  performance  of  the  sub-band  modulating 
signals. 

Although  the  sonar  platform  and  the  scatterer  were  somewhat  similar  to  the 
ones  simulated  in  section  4.4,  they  also  exhibited  some  key  differences,  which  inevitably 
influenced  the  comparison  between  the  simulated  data  and  the  in-water  data.  The  trans¬ 
mitter/receiver  on  the  sonar  platform  was  a  thirty-seven  element  array  of  transducers 
instead  of  a  single  omni-directional  transducer,  and  the  receive  beam  of  the  array  was 
steered  to  track  the  scatterer  position.  This  steering  was  done  both  physically  during 
the  data  collection  and  with  beamforming  techniques  dtuing  the  processing,  to  ensure 
that  the  scatterer  fell  within  the  mainlobe  of  the  receive  beam.  This  provided  the  high¬ 
est  SIR.  and  chance  for  detection.  The  weighting  and  filtering  for  the  received  array 
data  was  completed  in  Matlab.  While  the  steering  of  the  receive  beam  increased  the 
probability  of  detection,  the  geometry  of  the  experiment  resulted  in  exceptionally  small 
Doppler  values.  The  sonar  platform  was  stationary  over  all  of  the  interrogations,  while 
the  scatterer  always  produced  a  small  Doppler  shift  in  the  backscattered  signal.  A  carrier 
without  a  frequency  shift  was  used  in  the  demodulation  process,  because  the  scatterer 
induced  Doppler  was  close  to  zero  for  most  of  the  interrogations,  and  because  the  SIR^ 
was  so  low,  that  the  scatterer  return  could  not  be  discerned  from  the  rest  of  the  time  se¬ 
ries  for  frequency  matching.  Attempts  at  frequency  matching  were  made,  but  ultimately 
the  first  method  ended  up  being  more  effective. 

The  distance  between  the  scatterer  and  sonar  platform  was  much  larger  in  the 
in- water  portion  of  the  research  than  in  the  simulations,  which  resulted  in  smaller  signal- 
to-interference  ratio  measmrements  at  the  receiver.  The  scatterer  was  also  at  a  greater 
depth  in  the  ocean  than  the  sonar  platform,  so  the  temperature  and  salinity  strata  in 
the  ocean  had  a  greater  effect  on  the  propagating  signal,  which  was  not  taken  in  account 
in  the  simulations.  Several  attempts  were  made  to  better  model  the  in-water  data  in  the 
synthetic  data  generation  process  by  implementing  more  complex  simulations.  This  did 
not  significantly  aid  the  analysis  of  the  in-water  data,  or  provide  additional  insight  into 
the  performance  of  the  sub-band  modulated  signals  in  the  in-water  scenarios,  so  further 
simulations  were  suspended. 

Despite  the  differences  between  the  in-water  experiment  and  the  simulations,  the 
results  verified  that  sub-band  modulation  can  be  accomplished  in  realistic  situations. 
The  analysis  of  the  in-water  data  serves  as  a  confirmation  of  the  results  obtained  from 
the  simulations.  The  results  of  this  comparison  are  presented  and  analyzed  in  chapter 
5.2. 
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Chapter  5 

Results  and  Analysis 

5.1  Simulation  Analysis 

An  analysis  of  the  detection  performance  of  sub-band  modulation  is  made  through 
the  comparison  of  matched  filter  responses  for  different  signal  combinations  and  scenar¬ 
ios.  First,  a  few  individual  matched  filter  responses  for  sub-band  modulated  and  equiv¬ 
alent  unmodulated  signals  are  compared  in  order  to  analyze  the  sub-band  modulation 
process’s  effect  on  matched  filter  peak  shape,  location  and  detection  resolution.  Next, 
the  average  matched  filter  peak  responses  versus  signal-to-interference  ratio  at  the  re¬ 
ceiver  for  sub-band  modulated  and  unmodulated  signals  axe  investigated,  and  the  effects 
of  interference  on  scatterer  detection  for  each  process  are  compared. 

The  matched  filter  responses  for  several  of  the  sub-band  modulated  signals  were 
compared  to  the  matched  filter  responses  of  the  equivalent  unmodulated  signals  to  ensure, 
that  the  signals  were  not  distorted  by  the  sub-band  modulation  process.  The  goal  of  this 
analysis  is  to  verify  that  the  same  delay  and  Doppler  information  can  be  gathered  with  the 
sub-band  modulated  signals,  that  is  gathered  with  the  equivalent  unmodulated  signals. 
Figure  5.1  compares  the  matched  filter  response  of  the  sub-band  modulated  and  demod¬ 
ulated  BPSK  and  CW  signals  from  simulation  set  [®] 

and  s  ,  „ 
uh  2 

scatterer  at  an  SIR.  of  20dB.  The  carrier  used  with  the  discriminator  matched  the  re- 

t 

ceived  carrier  frequency  in  this  case.  Both  the  sub-band  modulated  BPSK  and  CW 
signals  estimated  identical  scatterer  delay  and  Doppler  parameters  to  the  equivalent 
unmodulated  signals.  Furthermore,  the  delay  and  Doppler  parameters  estimated  with 
the  matched  filters  accurately  represented  the  scatterer  parameters  set  up  in  this  simula¬ 
tion.  The  overall  phase-plane  shapes  of  the  matched  filter  response  peaks  are  also  almost 
identical,  which  indicates  that  the  sub-band  modulation  did  not  significantly  distort  the 
sonar  signals.  The  particular  replica  maps  used  in  this  portion  of  the  thesis  were  chosen, 
because  they  provided  a  large  enough  delay-Doppler  range  to  observe  the  matched  filter 
response  in  a  broader  section  of  the  phase-plane  than  just  the  area  immediate  to  the 
peak  response.  Although  this  did  slightly  compromise  the  detection  resolution,  it  was 
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[b]  BPSK  (unmod.)  max  =  -0.038904dB 
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[c]  CW  (sub.)  max  =  -0.8021 3dB 
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Fig.  5.1.  Comparison  between  Matched  Filter  Responses  of  Sub-band  Modulated  and 
Unmodulated  BPSK  and  CW  Signals  at  20dB  SIR^  (stationary  scatterer  and  matched 
carrier):  [a]  Sub-band  Modulated  BPSK  in  ?  [b]  Equivalent  Un¬ 
modulated  BPSK,  [c]  ^[t)  Sub-band  Modulated  CW  in  [d]  ^{t)  Equiv¬ 

alent  Unmodulated  CW 


65 


done  to  ensure  that  no  abnormalities  were  generated  in  the  matched  filter  response  by 
the  sub-band  modulation  process.  The  only  difference  in  the  matched  filter  responses  of 
the  sub-band  modulated  signals  of  simulation  set  a-nd  the  responses  of  the  corre¬ 
sponding  unmodulated  signals  is  a  decrease  in  the  peak  response  level  of  the  sub-band 
modulated  signals  by  approximately  IdB.  This  discrepancy  represents  a  drop  in  corre¬ 
lation  from  ~  1  to  ~  0.9  between  the  received  signal  and  the  replicas.  This  indicates, 
that  the  sonar  signals  are  slightly  distorted  by  the  sub-band  modulation  process,  but  the 
amount  of  distortion  is  minimal.  It  is  also  important  to  take  into  account,  that  the  signal 
energy  was  split  between  the  sub-band  modulating  signals.  While  both  the  modulated 
and  unmodulated  signals  had  approximately  the  same  SIR^  levels  when  received,  the 
signal  energy  is  split  between  the  individual  sub-band  signal,  once  they  axe  demodulated. 
Therefore,  the  signal-to-interference  ratio  is  lowered  for  each  of  the  demodulated  signals. 
Furthermore,  it  will  be  shown  later  in  this  section,  that  detection  performance  decreases 
as  the  signal-to-interference  ratio  decreases.  When  this  is  taken  into  consideration,  it  is 
not  unreasonable  to  expect  a  decrease  in  matched  filter  level  of  the  sub-band  modulated 
signals,  as  compared  to  the  unmodulated  responses  at  the  same  SIR^. 

The  matched  filter  responses  of  these  same  signals  are  also  compared  for  a  closing 
scatterer  at  3  (m/sec.)  and  an  opening  scatterer  at  -3  (m/sec.)  (figures  5.2a-d  and  5.3a- 
d  respectively).  The  same  carrier  used  above  was  implemented  in  these  demodulation 
processes,  which  resulted  in  a  carrier  mismatch.  Although  both  sets  of  figures  indicate 
that  sub-band  modulated  signals  estimate  the  same  delay  and  Doppler  parameters  of  a 
scatterer  as  the  equivalent  unmodulated  sonar  signals,  the  differences  between  peak  lev¬ 
els  of  the  matched  filter  responses  for  the  sub-band  modulated  signals  and  the  equivalent 
unmodulated  signals  are  much  greater  than  the  differences  seen  in  figures  5.1a-d  for  the 
stationary  scatterer.  The  matched  filter  peaks  for  the  moving  scatterers  were  approxi¬ 
mately  6dB  lower  thaxi  the  peaks  for  the  stationary  scatterer  (or  a  drop  in  correlation  to 
~  0.5).  This  decrease  in  peak  matched  filter  level  is  a  result  of  the  frequency  mismatch 
between  the  carrier  of  the  received  echo  signal  from  a  moving  scatterer  and  the  carrier 
implemented  in  the  demodulation  process,  which  further  distorts  the  received  signals. 

A  series  of  simulations  were  completed  to  characterize  the  performance  of  sub¬ 
band  modulation,  when  a  mismatched  carrier  is  implemented  in  the  demodulation  pro¬ 
cess  due  to  a  moving  scatterer.  A  series  of  simulations  were  run  using  the  experimental  set 
up  discussed  in  section  4.4.  The  scatterer  velocity  ranged  from  9(m/sec.)  to  -9(m/sec.), 
and  all  of  the  data  used  was  taken  at  7dB  SIR^  (±1.5dB).  All  of  the  simulations  were 
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Fig.  5.2.  Compaxison  between  Matched  Filter  Responses  of  Sub-band  Modulated  and 
Unmodulated  BPSK  and  CW  Signals  at  20dB  SIR^  (closing  scatterer  at  3  (m/sec.) 
and  mismatched  caxrier):  [a]  ^{t)  Sub-band  Modulated  BPSK  in  j  [b]  ^uh  1^^^ 

Equivalent  Unmodulated  BPSK,  [c]  Sub-band  Modulated  CW  in  si^d  [d] 

5^^2(^)  Equivalent  Unmodulated  CW 
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[c]  CW  (sub.)  max  =  -5.9293dB 
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Fig.  5.3.  Compaxison  between  Matched  Filter  Responses  of  Sub-band  Modulated  and 
Unmodulated  BPSK  and  CW  Signals  at  20dB  SIR.  (opening  scatterer  at  -3  (m/sec.) 
and  mismatched  carrier):  [a]  ^(t)  Sub-band  Modulated  BPSK  in  >  [b]  ^(t) 

Equivalent  Unmodulated  BPSK,  [c]  Sub-band  Modulated  CW  in  and  [d] 

^uh  2^^^  Equivalent  Unmodulated  CW 
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demodulated  using  the  unshifted  carrier  frequency  (w^),  which  produced  matched  fil¬ 
ter  responses  for  a  wide  range  of  mismatched  carrier  frequencies.  Three  simulations 
resulting  in  matched  filter  detections  were  computed  at  each  of  the  different  scatterer 
velocities,  and  each  set  of  three  peak  matched  filter  levels  were  averaged.  This  was  done 
to  ensure  that  the  typical  peak  level  of  the  matched  filter  output  and  not  an  atypical 
response  was  analyzed.  The  averaged  peak  responses  were  then  plotted  versus  the  differ¬ 
ence  between  the  transmitted  and  received  carrier  frequencies.  Figure  5.4  shows  a  6dB 
difference  between  the  peak  matched  filter  response  when  the  received  and  demodulated 
carrier  frequencies  are  identical  and  when  they  are  not.  This  equates  to  a  drop  in  the 
correlation  between  the  received  signal  and  the  replicas  from  '^.9  to  ~.5.  The  figure  also 
shows  that  the  performance  of  the  sub-band  signals  demodulated  with  a  mismatched 
carrier  does  not  vary  significantly  as  the  frequency  difference  increases.  Additional  sim¬ 
ulations  were  conducted  over  a  range  of  small  scatterer  velocities  near  0  (m/sec.)  to 
establish  a  frequency  mismatch  for  the  demodulation  process.  It  was  found,  that  a  range 
of  frequencies  1.5Hz  above  and  below  the  carrier  (or  Aw  =  ±0.2387  radians)  produced 
detection  peaks  equivalent  to  the  frequency  matched  response.  This  also  corresponded 
to  the  narrowband  Doppler  resolution  of  the  signal  (2.35).  Therefore,  carrier  frequency 
mismatches  falling  within  this  tolerance  can  still  be  considered  matched.  Another  series 
of  simulations  was  demodulated  with  a  carrier,  that  matched  the  received  carrier  fre¬ 
quency  from  a  scatterer  closing  at  3(m/sec.).  The  results  from  these  simulations  were 
identical  to  those  observed  in  figure  5.4.  A  prominent  peak  was  observed  at  the  matched 
frequency  and  the  performance  at  the  mismatched  frequencies  were  approximately  6dB 
lower.  Therefore,  sub-band  demodulation  yields  more  desirable  results  when  the  received 
signal  is  demodulated  with  a  matching  caxrier  (within  the  tolerance),  but  it  does  yield 
detection  results  for  mismatched  carriers  and  that  performance  was  relatively  constant 
over  the  range  of  scatterer  velocities  investigated.  These  results  are  also  displayed  for 
the  BPSK  and  CW  in  ^’^(t)  figure  5.5,  where  the  performance  of  the  sub-band  sig¬ 
nals  demodulated  with  a  mismatched  carrier,  the  sub-band  signals  demodulated  with  a 
matched  carrier  and  the  equivalent  unmodulated  signals  are  compared. 

Similar  detection  performance  was  observed  when  comparing  the  signals  of  '^^(t), 
'’PJt)  and  '0^(t)  with  corresponding  unmodulated  signals.  Examples  of  this  performance 
for  the  LFM  and  CW  signals  of  V’j,(t)  against  a  stationary  scatterer  are  seen  in  figures 
5.6a-d.  All  of  the  sub-band  modulated  signals  investigated  yielded  delay  and  Doppler 
information  consistent  with  the  scatterer  information  estimated  with  the  corresponding 
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Carrier  Mismatch  Results  for  BPSK and  CW  Signals  of  \|;  (t)  at  7dB  SIR. 

di  I 


Fig.  5.4.  Carrier  Mismatch  Results  for  BPSK  and  CW  Signals:  Magnitude  of  Matched 
Filter  Peaks  vs.  Frequency  Difference  between  the  Received  Carrier  and  the  Carrier 
Implemented  in  the  Demodulation  Process 
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[a]  BPSK  (sub.  mismatch)  max  =  -6.493ciB 
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[e]  BPSK  (unmod.)  max  =  -0.741 53dB 


[f]  CW  (unmod.)  max  =  -0.77833dB 
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Fig.  5.5.  Comparison  between  Matched  Filter  Responses  of  BPSK  and  CW  Signals 
at  7dB  SIR.:  Sub-band  Modulated  with  Matched  Carrier,  Sub-band  Modulated  with 

I 

Mismatched  Carrier  and  Unmodulated  (closing  scatterer  at  3(m/sec.)  ):  [a] 
Sub-band  Modulated  BPSK  with  Mismatched  Carrier  (Large  Mismatch)  ,  [b] 
Sub-band  Modulated  CW  with  Mismatched  Carrier,  [c]  Sub-band  Modulated  BPSK  with 
Matched  Carrier  (within  Tolerance)  ,  [d]  Sub-band  Modulated  CW  with  Matched  Carrier, 
[e]  j(i)  Equivalent  Unmodulated  BPSK,  and  [f]  ^{t)  Equivalent  Unmodulated 
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[a]  LFM  (sub.)  max  =  -0.36495dB 
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Fig.  5.6.  Comparison  between  Matched  Filter  Responses  of  Sub-band  Modulated  and 

Unmodulated  LFM  and  CW  Signals  at  20dB  SIR,  (stationary  scatterer  and  matched 

% 

carrier):  [a]  ^{t)  Sub-band  Modulated  LFM  in  ,  [b]  ^{t)  Equivalent  Unmod¬ 
ulated  LFM,  [c]  5^^  ^{t)  Sub-band  Modulated  CW  in  and  [d]  ^{t)  Equivalent 

Unmodulated  CW 
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unmodulated  signals.  They  also  exhibited  the  same  decreased  matched  filter  perfor¬ 
mance  as  the  signals  of  ip  {t),  when  demodulated  with  a  mismatched  carrier.  The  only 
significant  difference  observed  was  a  slight  decrease  in  the  peak  matched  filter  levels  for 
the  CW  signals  in  simulation  sets  and  by  ~ldB  in  both  the  matched  and 
mismatched  carriers  cases.  A  decrease  was  to  be  expected,  since  the  signal  energy  is  split 
between  four  sonar  signals  in  these  sets,  instead  of  only  two  as  in  sets  and 

Therefore,  it  can  be  concluded,  that  the  peak  matched  filter  response  for  each  signal 
decreases  according  to  the  number  of  sub-band  signals  modulating  a  given  carrier,  but 
their  responses  are  otherwise  undistorted.  However,  all  of  the  signals  used  in  each  of  the 
simulations  sets  behaved  similarly  for  the  most  part,  which  illustrates  that  the  sub-band 
modulation  detection  performance  is  independent  of  signal  type,  but  is  dependent  on 
number  of  modulating  signals  to  a  certain  degree. 

Although  demodulating  received  sub-band  signals  with  matching  carriers  yields 
stronger  matched  filter  responses  than  those  of  signals  demodulated  with  mismatched 
carriers  (figure  5.5),  matching  the  received  carrier  frequency  with  an  ideal  discrimina¬ 
tor  presents  other  problems.  There  is  no  a  priori  knowledge  of  the  carrier  frequency, 
so  it  needs  to  be  measured  in  order  to  be  matched.  Frequency  matching  (or  frequency 
tracking)  of  the  received  signal  requires  an  estimation  of  the  Doppler-shifted  caiider 
frequency,  which  can  be  very  time  consuming  and  computationally  intensive.  This  is 
especially  difficult,  when  the  SIR.  is  low.  The  effects  of  matched  versus  mismatched 
carrier  demodulation  on  detection  performance  coupled  with  frequency  matching  diffi¬ 
culties  prompted  an  investigation  into  the  implementation  of  the  demodulation  process 
as  part  of  the  hypothesis  testing.  The  received  signal  was  demodulated  with  an  array  of 
time  dilated  carriers,  as  discussed  in  section  3.3.1,  and  the  resulting  array  of  demodulated 
signals  were  spectrally  separated.  The  arrays  of  signals  were  matched  filtered,  and  an 
improvement  in  detection  performance  was  observed  (figure  5.7).  Only  a  small  number  of 
simulations  were  completed  using  this  method  of  demodulation  and  hypothesis  testing, 
but  they  did  yield  positive  results.  Since  the  coupled  demodulation,  hypothesis  testing 
process  implements  frequency  matching,  it  is  not  surprising,  that  resulting  matched  filter 
responses  had  maximum  values  around  the  same  level  as  the  matched  carrier  responses 
seen  earlier  in  this  analysis.  The  detection  performance  observed  in  these  experiment 
demonstrate  that  a  further  investigation  into  the  implementation  of  demodulation  as 
part  of  hypothesis  testing  would  be  a  good  topic  for  future  research.  This  performance 
also  indicates,  that  it  may  be  worthwhile  to  investigate  demodulators  with  frequency 
tracking,  such  as  the  phase  locked  loop  discussed  in  section  3.3.1. 


73 


[a]  Demodulation  Coupled  with  Hypothesis  Testing  (BPSK)  max  =  -1-2729dB 
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[b]  Demodulation  Coupled  with  Hypothesis  Testing  (CW)  max  =  -1.736dB 
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Fig.  5.7.  Matched  Filter  Responses  of  Sub-band  Modulated  BPSK  and  CW  signal 
Combination  sit  7dB  SIR^,  when  Demodulation  is  Coupled  with  Hypothesis 

Testing:  [a]  Sub-band  Modulated  BPSK  and  [b]  Sub-band  Modulated 

CW 
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The  effects  of  noise  and  interference  on  the  detection  performance  of  the  sub-band 
modulating  signals  were  investigated  through  the  analysis  of  the  resulting  matched  filter 
responses  at  different  SIR^  levels.  A  series  of  simulations  were  conducted  over  a  range  of 
SIR.  values  from  20dB  to  -13dB  with  the  BPSK  and  CW  signals  in  'ip  {t) .  Ten  detections 

I  O' 

were  made  at  each  SIR.,  and  their  peak  matched  filter  response  levels  were  then  aver- 

aged.  These  average  peak  levels  for  the  matched  filters  were  then  plotted  with  respect 

to  SIR..  This  was  done  with  a  stationary  scatterer,  and  then  repeated  with  a  scatterer 

closing  at  3(m/sec.)  using  the  same  cajrier  for  demodulation  to  observe  the  effects  of  a 

mismatched  carrier  on  the  noise  performance  of  the  sub-band  signals  (figure  5.8).  Since 

the  scatterer  parameters  are  controlled  in  these  simulations,  it  was  straight  forward  to 

verify  the  accuracy  of  the  parameter  estimations  made  with  the  matched  filters.  This 

was  done,  so  only  accurate  detections  were  used  as  the  ten  peak  matched  filter  responses 

averaged  at  each  SIR.  investigated.  Therefore,  the  effects  of  interference  on  the  perfor- 

mance  of  the  sub-band  modulating  signals  could  be  assessed  independently,  instead  of 

coupled  with  the  probability  of  detection.  A  proper  investigation  of  the  probability  of 

detection  for  the  ideal  discriminator  would  be  more  meaningful,  when  compared  with 

the  probability  of  detection  of  several  other  demodulators  This  comparison  is  another 

topic  left  for  future  research. 

Figure  5.8  shows  a  decrease  in  the  peak  matched  filter  responses  as  the  signal- 
to-interference  ratio  at  the  receiver  decreases.  Since  the  curves  for  the  BPSK  and  CW 
signals  were  essentially  the  same  for  each  case,  they  were  combined  to  provide  a  more 
general  description  of  the  noise  performance  of  sub-band  modulation  independent  of 
signal  type.  The  curves  for  the  sub-band  modulated  (both  the  matched  and  mismatched 
carrier  cases)  and  equivalent  unmodulated  signals  followed  roughly  the  same  shape,  only 
at  different  levels.  The  relative  levels  of  the  three  curves  in  figure  5.8  corresponded  to 
the  peak  matched  filter  levels  observed  earlier  in  this  analysis.  The  average  response 
of  the  sub-band  signals  demodulated  with  a  matched  carrier  is  about  IdB  lower  than 
the  unmodulated  response  for  SIR^  values  above  OdB  (or  a  drop  in  correlation  of  ~ 
0.1).  Accordingly,  the  average  response  of  the  sub-band  signals  demodulated  with  a 
mismatched  carrier  (outside  of  the  tolerance)  is  about  6dB  lower  than  the  response  of 
the  matched  carrier  case  (or  7dB  below  the  unmodulated  response)  for  RIR^  values 
above  OdB.  This  corresponds  to  a  drop  in  correlation  of  ~  0.4  below  the  matched  carrier 
case  and  ~  0.5  below  the  unmodulated  case.  All  three  cases  were  relatively  indifferent 
to  interference  above  OdB  SIR^,  but  the  peak  matched  filter  responses  began  to  steadily 
decrease  as  the  SIR.  decreased  further.  This  closely  follows  the  demodulation  theory 


Average  of  the  Peak  Matched  Filter  Responses  of  Signals  In  over  240  Runs  with  Respect  to  SIR. 
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Fig.  5.8.  Average  of  the  PeaJc  Matched  Filter  Responses  of  Signals  in  over  240 

Simulations  with  Respect  to  SIR.  for  Three  Cases:  Sub-band  Signal  Demodulated  with 

% 

a  Matched  Carrier  (within  the  Tolerance),  Sub-band  Signal  Demodulated  with  a  Mis¬ 
matched  Carrier  (Large  Mismatch)  and  the  Equivalent  Unmodulated  Signal 
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presented  in  section  3.3.2.  Similar  results  were  obtained  when  other  simulation  sets 
and  scatterer  velocities  were  used  Therefore,  the  results  displayed  in  figure  5.8  are  a 
reasonable  representation  of  the  overall  noise  performance  of  sub-band  modulation. 

5.2  In- Water  Data  Analysis 

The  matched  filter  responses  of  the  in-water  data  were  analyzed  to  verify  the 
detection  performance  of  sub-band  modulation  in  a  real  world  scenario.  However,  this 
was  more  difficult  than  anticipated.  The  reverberation  levels  for  each  interrogation  were 
so  high,  that  the  scatterer  returns  were  not  discernible  in  most  of  the  time-series.  This 
resulted  in  very  low  SIR.  levels  at  the  receiver.  The  strongest  SIR.  level  detected  was 
approximately  OdB,  and  the  average  SIR^  values  were  below  -5dB  for  the  majority  of 
the  interrogations. 

Each  of  the  twenty-four  scatterer  returns  were  demodulated  and  matched  filtered 
against  the  corresponding  replica  maps.  The  peaks  in  matched  filter  responses  were  at  a 
relatively  low  level,  since  the  SIR^  was  so  low  for  each  of  the  interrogations.  The  peaks 
of  the  matched  filter  responses  for  the  BPSK  signals  were  noticeably  lower  than  those 
of  the  CW  signals,  even  though  the  matched  filter  peaks  for  the  sub-band  modulated 
CW  and  BPSK  signals  were  roughly  the  same  level  observed  in  the  simulations.  Fur¬ 
thermore,  the  level  of  the  BPSKs’  responses  were  low  enough,  that  significant  peaks  in 
the  matched  filter  responses  were  only  detected  in  nine  interrogations  out  of  twenty-four 
(and  only  four  of  those  detections  were  prominent  and  easily  discernible  in  the  phase 
plane),  whereas  significant  matched  filter  peaks  were  detected  in  sixteen  interrogations 
for  the  CW  signals.  Seven  of  these  CW  detections  were  made  without  corresponding 
BPSK  detections,  while  the  other  nine  were  from  the  same  interrogations  as  the  BPSK 
detections  (otherwise  known  as  joint  detections).  No  BPSK  detections  were  made  with¬ 
out  corresponding  CW  detections.  This  difference  in  detection  performance  with  the 
BPSK  signals  is  the  result  of  it’s  phase  sensitivity  combined  with  the  low  SIR^  levels  in 
these  interrogations.  Despite  this  difference,  characteristics  of  the  sub-band  modulating 
signals’  performance  can  still  be  analyzed  from  the  few  matched  filter  detections. 

Sixteen  detections  were  made  in  twenty-four  interrogations  with  the  sub-band 
modulated  CW  signal.  The  delay  and  Doppler  values  estimated  from  these  matched 
filter  responses  corresponded  with  the  scatterer  parameters  observed  during  the  data 
collection  for  the  majority  of  these  interrogations.  Nine  of  these  interrogations  were 
also  joint  detections  with  both  the  CW  and  BPSK  signals,  and  the  Doppler  estimations 
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from  both  of  the  sub-band  signals  were  equivalent.  The  delay  estimations  were  also 
approximately  equivalent,  except  that  the  BPSK’s  delay  estimate  was  at  a  much  higher 
resolution,  which  was  expected  (section  2.2).  Figure  5.9  shows  example  matched  filter 
responses  for  one  of  these  interrogations,  that  resulted  in  a  joint  detection.  Only  eight 
interrogations  out  of  twenty-four  did  not  result  in  scatterer  detections  with  either  the  sub¬ 
band  modulated  CW  signals  or  both  the  sub-band  modulated  CW  and  the  BPSK  signals. 
The  matched  filter  response  for  these  interrogations  were  predominantly  obscured  by 
reverberation. 

Despite  the  low  number  of  joint  detections,  the  in-water  sub-band  modulated 
signals  performed  reasonably  similar  to  those  in  simulations  with  low  5/i?/s,  that  were 
also  demodulated  with  mismatched  carriers.  This  was  especially  reassuring,  since  the 
SIR.  was  too  low  to  effectively  frequency  match  the  received  signals  in  the  in- water 
case.  Both  the  in-water  data  and  the  simulations  verified,  that  multiple  sonar  signals 
sub-band  modulating  a  carrier  could  be  used  to  simultaneously  estimate  a  scatterer’s 
parameters.  They  also  showed,  that  sub-band  modulated  signals  could  still  be  used  to 
detect  a  scatterer  even  at  low  SIR.  values  and  without  precisely  matching  the  received 
carrier  frequency.  Although  there  were  only  a  limited  number  of  in-water  interrogations 
available  for  analysis  in  this  thesis,  they  yielded  favorable  results,  even  though  they  were 
at  very  low  SIR.  values  against  a  scatterer  with  very  small  l.o.s.  velocities.  Therefore, 
the  further  investigation  of  sub-band  modulation  through  more  extensive  in- water  testing 
against  a  variety  of  scatterers  with  different  Dopplers  over  a  large  range  of  SIR^  values 
might  be  a  worthwhile  topic  for  future  research. 
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Fig.  5.9.  Matched  Filter  Responses  of  Sub-band  Modulated  BPSK  and  CW  Signal 
Combination  (■;/;  (t))  (In-Water):  [a]  ^  (t)  Sub-band  Modulated  BPSK  and  [b]  2(^) 

Sub-band  Modulated  CW 
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Chapter  6 


Summary  and  Conclusions 


This  thesis  presented  a  method  of  simultaneously  transmitting  multiple  sonar  sig¬ 
nals,  while  maximizing  transmission  time  and  avoiding  losses  in  total  transmit  energy 
caused  by  amplitude  modulation  in  a  peak  power  limited  system.  This  process  is  known 
as  sub-band  modulation,  and  it  is  the  linear  combination  of  multiple  sonar  signals  as 
the  modulating  phase  component  of  a  carrier  signal.  An  analysis  of  the  advantages  and 
disadvantages  of  different  sonar  waveforms  in  different  situations  illustrated  the  benefits 
of  transmitting  multiple  sonar  signals.  The  specific  benefits  of  different  sonar  signal 
combinations  were  also  discussed.  An  extensive  background  of  sonar  theory  and  echo 
processing  explained  the  advantages  of  simultaneously  transmitting  multiple  sonar  sig¬ 
nals  to  characterize  a  given  underwater  event  and  the  need  to  limit  transmission  time  in 
active  sonar.  This  background  also  gave  insight  into  the  methods  used  to  estimate  delay 
and  Doppler  parameters  of  a  scatterer  from  echo  signals  backscattered  off  of  it.  Then, 
different  types  of  modulation  were  investigated  and  compared  to  show  the  disadvantages 
of  amplitude  modulation  in  a  peak  power  limited  system  and  to  further  justify  the  use 
of  sub-band  modulation. 

Next,  the  theory  behind  sub-band  modulation  was  presented.  The  factors  govern¬ 
ing  the  modulation  were  discussed,  and  several  examples  were  examined,  that  illustrated 
effects  of  these  factors  on  transmit  energy,  modulation  bandwidth  and  sideband  loca¬ 
tions  of  the  modulating  sub-bands.  This  was  followed  by  an  analysis  of  the  demodulation 
process  used  to  extract  the  individual  sonar  signals  from  the  backscattered  echo  signals. 
Special  attention  was  spent  on  the  demodulation  difficulties  due  to  the  scatterer  pa¬ 
rameter  effects  on  the  echo  signal.  Most  of  these  difficulties  stemmed  from  a  frequency 
mismatch  between  the  carrier  of  the  received  echo  signal  and  the  carrier  implemented 
in  the  demodulation  process.  Methods  of  resolving  these  problems  through  frequency 
matching  and  coupling  the  demodulation  process  with  hypothesis  testing  were  then  dis¬ 
cussed  and  analyzed.  These  methods  were  found  to  successfully  counteract  the  carrier 
mismatch.  This  discussion  was  followed  by  an  analysis  of  the  effects  of  interference  on 
the  demodulation  of  a  sub-band  modulating  signal.  In  this  section,  it  was  found,  that 
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the  performance  of  the  sub-band  demodulation  remains  relatively  unaffected  by  inter¬ 
ference,  when  the  signal-to-interference  ratio  is  high.  When  the  SIR^  drops  low  enough, 
the  demodulator  performance  rapidly  degrades. 

The  practicality  of  sub-band  modulation  was  discussed  in  the  experiment  section 
of  this  thesis.  First,  waveform  design  considerations  were  discussed,  and  methods  of 
selecting  the  modulation  index  and  the  frequency  spacing  of  sub-band  waveforms,  while 
taking  into  consideration  transducer  bandwidth  and  potential  scatterer  Doppler,  were 
presented.  Second,  the  different  sub-band  signal  combinations  implemented  in  the  ex¬ 
perimental  portion  of  this  thesis  were  presented,  and  the  motivation  behind  these  signal 
combinations  were  discussed.  An  analysis  of  the  modulation  index  used  in  this  research 
followed  this  discussion.  Third,  the  Matlab  implementation  of  the  sub-band  modulation 
and  demodulation  processes  were  summarized,  along  with  the  echo  processing.  Finally, 
an  overview  of  the  different  simulated  and  in- water  experiments  completed  in  this  thesis 
and  their  goals  were  discussed.  These  experiments  investigated  the  overall  detection 
performance  of  sub-band  modulating  signals  and  were  conducted  to  examine  the  effects 
of  modulating  signal  type,  signal-to-interference  ratio  and  demodulation  conditions  on 
detection. 

Finally,  the  results  of  the  simulated  and  in-water  experiments  were  discussed  and 
analyzed.  For  the  most  part,  the  detection  performance  of  the  sub-band  modulating 
signals  was  investigated  in  the  simulations,  while  the  in-water  experiment  was  used  to 
confirm  this  performance.  The  simulations  indicated,  that  the  signal  type  or  combination 
does  not  adversely  affect  the  detection  performance  of  the  sub-band  modulating  signals, 
while  the  number  of  sub-band  modulating  signals  did  influence  detection.  Scatterer 
range  and  velocity  estimations  made  from  the  matched  filter  responses  of  the  sub-band 
modulating  signals  matched  the  scatterer  parameters  in  the  simulations,  and  the  shapes 
of  the  matched  filter  responses  were  not  noticeably  different  from  the  responses  of  equiv¬ 
alent  unmodulated  sonar  signals.  However,  the  peak  magnitudes  of  the  matched  filter 
responses  of  the  sub-band  modulating  signals  were  slightly  lower  than  the  responses  of 
the  corresponding  unmodulated  signals,  because  of  a  small  drop  in  correlation  between 
the  demodulated  echo  signal  and  the  bank  of  replicas  caused  by  the  sub-band  modulation 
process.  It  was  also  determined,  that  the  detection  performance  of  the  sub-band  signals 
were  negatively  was  affected  by  carrier  mismatch.  The  demodulation  process  used  in  this 
thesis  works  best,  when  the  demodulator  can  match  the  carrier  frequency  of  the  received 
signal.  This  performance  was  confirmed  by  the  results  from  the  simulations  involving 
the  coupled  demodulation,  hypothesis  testing  process.  The  performance  of  the  sub-band 
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modulating  signals  demodulated  with  a  matched  carrier,  the  sub-band  signals  demodu¬ 
lated  with  a  mismatched  carrier  and  the  equivalent  unmodulated  signals  were  evaluated 
and  compared  over  a  range  of  signal-to-interference  ratios  at  the  receiver.  Although 
there  was  a  difference  between  the  matched  filter  peak  magnitudes  in  the  three  cases, 
all  of  the  cases  showed  similar  decreases  in  performance  as  the  signal-to-interference  ra¬ 
tio  decreased.  A  small  set  of  in-water  data  with  a  sub-band  modulated  signal  was  also 
analyzed  to  verify  the  performance  seen  in  the  simulations.  The  results  of  this  analysis 
were  similar  to  the  performance  of  the  simulated  data  in  similar  circumstances,  which 
consisted  of  small  signal-to-interference  ratios  and  mismatched  carriers. 

Possible  future  research  topics  in  sub-band  modulation  could  included  a  further 
investigation  of  different  modulation  indices  and  the  resulting  modulation  bandwidths. 
Another  possibility  for  futmre  research  would  be  to  investigate  different  frequency  track¬ 
ing  demodulators,  such  as  the  phase  locked  loop,  and  compare  them  with  the  ideal 
discriminator.  Further  evaluation  of  the  coupled  demodulation,  hypothesis  testing  pro¬ 
cess  is  also  warranted.  Finally,  more  experiments  (both  simulated  and  in-water)  should 
be  conducted  using  more  complex  sceneries  and  a  variety  of  different  scatterers  resulting 
in  a  wide  range  of  Doppler  values  to  more  thoroughly  examine  the  detection  performance 
of  sub-band  modulating  signals. 
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